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ABSTRACT 


Various  authors  have  carried  out  studies  in  the 
general  area  of  the  measure  of  a  random  set*  W.L*  Stevens 
considered  arcs  of  length  D  (0  <  D  <  l)  at  random  points 
on  the  circumference  of  a  circle  of  length  unity  and 
derived  the  probability  that  the  circumference  is 
completely  covered*  A  somewhat  similar  but  non-circular 
problem  was  dealt  with  by  D.F.  Votaw  Jr*  He  derived  the 
probability  function  of  the  measure  of  a  random  linear  set. 
C*  Domb  also  considered  a  problem  of  random  intervals  on 
the  line*  H*E.  Robbins  derived  general  formulae  for  the 
moments  of  the  measure  of  any  random  set*  The  results  of 
these  authors  and  others  are  reviewed  in  the  first  two 
chapters  of  this  thesis* 

The  third  chapter  is  devoted  to  the  derivation 
of  the  probability  function  of  the  portion  of  the 
circumference  of  a  circle  of  unit  length  which  is  covered 
by  arcs,  of  length  D  (0  <  D  <  1)  ,  at  random  points  on 
the  circumference*  The  first  two  moments  of  the  coverage 
are  also  derived* 
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CHAPTER  I 


INTRODUCTION 

Consider  n +  1  points  z^ (i =  0 , 1 , . . . ,n) 
selected  independently  and  at  random  from  the  interval 
(0,1)  ,  the  distribution  of  any  z.  being  the  uniform 

l 

distribution  with  distribution  function 


(1.1) 


F ( z )  =  z  (0  <  z  <  1)  . 


Let  the  n+1  points  be  ordered  in  ascending  order  of 
magnitude  as  x^(i=  0,1,  .  .  .n)  .  The  totality  of  all  cases 
in  which  two  of  the  are  equal  has  zero  probability,  and 

such  cases  can  be  excluded  without  affecting  the  problem. 
Associate  with  each  x^  an  interval  I  where 


(1*2)  I 


(x . ,x.+D ) 

^  l’  1  7 

k  (xi,l)  U  (O.x^D-l) 


if  x.+D  <1  (0  <  D  <  1  ;i=  0,1,. 

if  x  .+D  >1 
1 


Let  X  denote  the  random  set  which  is  the  point  set  sum  of  the 
n+1  intervals  {i}  and  let  |i(X)  be  its  measure.  The  range 
of  li(X)  is  D  <  |i(X)  <  m  ,  where  m  denotes  the  minimum  of 
1  and  (n  +  l)  D  .  Using  the  methods  of  D.F.  Votaw  Jr.  (1946), we 
show,  in  the  third  chapter,  that  the  probability  function 


. 


CM 

■ 
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f^(x)  of  the  random  variable  |i(X)  is  given  by 

(1.3)  f(x)=n'.  £  £  (-1)  _  (n+1 

j=l  r=0  (n-j ) !  ( j-1 ) !  3 

where  qD  <  x  <  (q+1 )D, (q=l , 2 ,M;M = minimum  (n,[l/D]))  , 
x  <  m  ,  and  where  the  symbol  [x]  represents  the  greatest 
integer  less  than  x  .  When  m  =  (n+l)D  ,  the  distribution 
function,  Fn(x)  >  ° f  |l(X)  has  a  discontinuous  jump  or  saltus 
point  corresponding  to  the  finite  probability  that  the  unit 
interval  will  consist  of  exactly  n+1  non-overlapping  intervals 
{i}  •  At  any  such  discontinuity  of  F  (x)  the  function  fn(x)  t 
the  probability  function  of  |i(X)  ,  is  not  defined  in  the  pure 
mathematical  sense. 


)  (n-J+1)  (l-.x)J  ‘1[x-D(  j+r)  ]n_:j , 
r 


The  expected  value  of  |i(X)  is  found  to  be 


(1.4) 


E(H(X))  =  1  -  ( 1-D) 


n+1 


while  its  variance  is  given  by 


"7T  (2(1-D)n+2  +n  (l-2U)n+2)  -  (1-D)2n+2 
n+z 


(1.5)  Var  (m(X))= 


for  0  <  D  <  i  , 


-75-  (1-D)n+2  -  (1-D)2n+2  for  \  <  D  <  1  . 
n+Z  * 


To  our  knowledge,  the  majority  of  the  research  which 
one  might  classify  under  the  general  heading  of  the  measure  of 


- 
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a  random  set  was  carried  out  during  the  decade  following  a 
1939  paper  by  W.L.  Stevens*  In  the  second  chapter  of  this 
thesis ,  we  present  a  review  of  some  of  the  work  of  various 
authors  who  have  conducted  studies  in  this  area .  The 
remainder  of  this  chapter  is  devoted  to  presenting,  in 
chronological  order  beginning  with  Stevens*  paper,  the  results 
obtained  by  these  authors* 

W.L.  Stevens  (1939)  solved  the  following  problem: 
consider  n+1  arcs  of  length  D  (0  <  D  <  l)  at  random  points 
on  a  circle  of  unit  circumference.  Choose  the  end  of  one  arc 
arbitrarily  as  the  origin,  and  measure  distances  anti¬ 
clockwise  around  the  circle,  the  location  of  each  arc  being 
represented  by  its  first  point,  i.e.  the  point  at  the 
clockwise  end.  What  is  the  probability  that  every  point  of 
the  circle  is  included  in  at  least  one  of  the  arcs,  i.e* 
what  is  the  probability  that  the  circle  is  completely  covered  ? 
It  is  clear  that  this  representation  is  equivalent  to  the 
representation  of  the  linear  problem  described  previously. 
Stevens  proved  that  the  required  probability  is  given  by 

(1.6)  1  -  (“^1)(l-D)n  +  (n21)(l-2D)n  -  ...  ±  (“^Jd-kD)"  , 

where  k  is  the  greatest  integer  less  than  1/D  .  Stevens 
also  calculated  the  frequency  distribution  of  the  number  of 


mb  ©  >  o  brt«  e  if  o  .  z'\±o  dxii  i  o  i  lu  ixo  &  no 

tad  3T  ,  Ij  >  *•  o  -O  'I  >  I  rl  "  ,  BiMxo  <•  j-  bnuo'tfc  •*••■  >o!b 


vI»^»lqnoo  i  i  dlsiio  add  darid  xd±Iidado*iq  •  ai  j  wiw 

. 

add  od  dnoJtaviupa  b!  aojtdsdnaaa*iq&*i  axrfd  darfd  *xaaIo  *i  dl 

. 


gaps  (there  is  said  to  be  a  gap  after  the  r  th  arc,  if  for 
a  distance  greater  than  D  beyond  the  first  point  of  this  arc 
there  is  no  other  arc)  • 

H.E.  Robbins  (1944)  considered  a  similar  but 
non-circular  problem  and  proved  the  following 

Theorem:  Let  X  be  a  random  Lebesgue  measurable  subset  of 

Euclidean  n-dimensional  space  E^  ,  with  measure  p(X)  .  For 
any  point  x  of  E^  let  p(x)  =  Pr(x  €  X)  »  Suppose 
furthermore  that  the  function 


(1.7) 


g(x,X) 


1 

0 


x  6  X  , 
x  $  X 


is  a  measurable  function  of  the  pair  (x,X)  .  Then  the 
expected  value  of  the  measure  of  X  will  be  given  by  the 
Lebesgue  integral  of  the  function  p(x)  over  E^  ,  i.e. 


(1*8)  E  (|i(X) )  =  J  p(x)  dp (x )  . 

En 

He  also  generalized  this  result  by  showing  that  the  m  th 
moment  of  p(X)  is  given  by 

(1.9)  E(pm(X))  =  J*  p(x1,x2,  .  .  .  ,xn)  dp  (xx,x2,  .  . .  ,xn)  , 


‘  T  •' 
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where 


*  •  • » 


x  )  denotes  Lebesgue  measure  in  E 


and  where 


(1.10)  p(xx,x25 


} x  )  =  P  (xn  6  X  and  x_  €  X 
n  r  1  2 


•  •  • 


and  x  6  X)  . 
n 


In  the  following  problem,  let  the  random  set  X  be 


defined  as  follows:  let  A^,a^(i =  1,2, . . . ,n)  and  6  be 
fixed  positive  numbers  such  that  a^  <  26  ;  let  R  denote  the 
n-dimensional  interval  consisting  of  all  points 

(x^,x2, • • • ,xn)  such  that  0  <  x^  <  ,  and  let  R*  denote 

the  larger  interval  for  which  -  6  <  x^  <  +  6  (and  also 

n 

its  measure  II  (A. +  26))  .  Let  a  fixed  number  N  of 


i=l 


intervals  with  sides  parallel  to  the  axes  be  chosen  indepen 
dently  with  the  density  function  for  the  center  of  each 


interval  constant  and  equal  to  l/R/  in  R*  .  The  set  X 


is  the  intersection  of  the  point  set  sum  of  the  N  intervals 
with  R  •  What  is  the  mean  and  variance  of  the  measure  of 
the  set  X  ?  In  a  second  paper,  H.E.  Robbins  (1945),  using 
the  main  result  of  his  1944  paper,  showed  that 


N 


(l.ii) 


n 


II  a . 

,  i 


where  r 


and  that 


..  in  £Hlr  iNufii 


. 


N 


„  ”n  al  / 

(1.12)  Var  ( |i(X) )  =  2p  f  ...  f  1- 


<an“An>  <al"Al> 


n 

2r  -  n  w . 
i=l  1 

r' 


n 


,  TT  ( w  •  +A .  -a  .  )  dw,  ...  aw 
l  i  i  9  1  n 


9  N  n  0  n  9  o'io  _  2N 

+ < 1  -  ?  >  { ±n=1  Ai  Ai -< Ai-ai>  )}-R  ^1-r) 


where  'the  symbol  <!  x  >  is  defined  by 


(1.13) 


<  X  > 


a 


if  x  >  0  , 
if  x  <  0  * 


In  this  paper,  Robbins  also  considered  the  more  general  case 
where  the  number  N  of  n-dimensional  intervals  is  not  consid¬ 
ered  as  being  fixed  but  is  taken  as  a  random  variable  capable 
of  assuming  the  values  0,1,2,...  with  respective  probabilities 
Po>Pi,P2^***  j  and  with  generating  function 


CO 


(1.14) 


<p(t)  =  £  PNt 


N 


N=0 


He  showed  that 

E(p(X))  =  R  {l -  <P  (1  - 


3 


■ 


(M.I) 


and  that 


n 


(1.16)  Var  (n(X))  =  2  J  ...  J1  «P [  1 

<*~n  n 


a-~A->  <al~Al> 


n 

2r  -  II  w. 
i=l 1 

R 


n 


II  ( w .  +A .  -a  .  )  dw,  *  •  *  dw 
l  i  l  1  n 


n  2  n 


{  n  A  .  -  n  (A?  -  <  A.  -  a.  >*)  \  V  (l-§£  )~R2cp2(l  ) 

^1*1  1  i=i  i  i  i  '  J  R  K 


Robbins  also  derived  formulae  for  the  corresponding  problem 
for  circles  in  the  plane. 

J.  Bronowski  and  J.  Neyman  (1945),  using  methods 
quite  different  from  those  of  Robbins,  considered  a  two- 
dimensional  problem  and  found  the  mean  and  variance  of  the 
measure  of  the  random  set  X  ,  where  X  was  taken  as  the 
point  set  sum  of  rectangles  in  the  plane  (with  sides  parallel 
to  the  axes).  Their  results  agreed  with  the  results  of 
Robbins  (1945)  for  n=2  (see  (1.11),  (1.12),  (1.15),  and 
(1.16))  * 

In  a  later  paper,  L.A.  Santalo  (1947),  using  the 
methods  of  Robbins,  found  the  mean  and  variance  of  the 
measure  of  the  random  set  X  consisting  of  the  point  set 
sum  of  N  rectangles  in  the  plane  but  having  variable 
position.  A  similar  problem  was  solved  for  n  dimensions 
with  the  N  rectangles  replaced  by  spheres  • 


1 


.  Jem  aniaxj  nw rcaH  islawomna  #L 


. 
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Returning  to  the  one-dimensional  case. 


D.F.  Votaw  Jr*  (1946)  considered  a  non-circular  problem  and 
derived  the  probability  function  of  the  random  variable  y 
defined  as  follows:  consider  a  random  sample  z^(i=l,2,*..,n) 
of  n  values  of  a  one-dimensional  random  variable  Z 
having  distribution  function  F(z)  •  Arrange  the  n  sample 
values  in  increasing  order  of  magnitude  as  x^(i=  1,2,  •  •  •  ,n) 
and  consider  the  set  consisting  of  the  intervals 


(1.17) 


Denote  its  measure  by  p. ( X)  *  Let  y  =  ji(X)  -  D  •  The 
range  of  y  is  0  <  y  <  m  ,  where  m  denotes  the  minimum 
of  1  and  (n-l) D  ♦  Votaw  showed  that  if  F(z)  is  given 
by  (1*1)  ,  the  probability  function  fn(y)  y  is  given 


by 


q 


(1.18)  f  (y)  =  n  2 

XI  • 


3=0  r=0 


qD  <  y  <  (q+l)D,  (q=0,l 


M;M  =  minimum 


y  <  m  • 


The  expected  value  of  y  is 


si  i « v  mobn*n  ^  ::  ox  mib-dno  a  ‘io  Pd  .  57 


X  «  ^•n®C 

' 


(1.19) 
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E(y)  =  (^1)  I1'  (1'D)n+1]  • 

C*  Domb  (1947)  considered  random  intervals  with  a 
constant  length  D  ;  whose  left  endpoints  are  random  events 
distributed  according  to  a  Poisson  distribution*  He  let 
W(x,P)  be  the  probability  that  the  measure  of  the  subset  of 
the  interval  (0,P)  which  is  covered  by  random  intervals 
does  not  exceed  x  •  He  then  obtained  the  p-multiplied 
Laplace  transform  of  W(x,(3)  as  the  solution  of  an  integral 
equation  of  the  convolution  type*  From  this  he  derived  many 
results*  In  particular,  he  determined  the  probability  that 
the  interval  (0,P)  is  completely  covered  (or  completely 
uncovered  )•  He  also  derived  an  expansion  of  W(x,P)  and 
the  moments  of  the  W(x,P)  distribution*  The  method  also 
yielded  analogous  results  for  the  problem  we  have  considered, 
namely  that  of  random  arcs  on  a  circle* 


•  oqvd’  noX&uj.oYaoo  arfcfr  no.  daxjpa 


'  ^cfrxIxdxK  <j  ori  fas  tacfrob  exi  1 


rlo  noianaqxe  xib  faovlieb  oala  aH 
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CHAPTER  II 


REVIEW  OF  SOME  RELATED  PROBLEMS 


As  was  mentioned  earlier,  this  chapter  is  intended 
to  furnish  the  reader  with  a  review  of  the  solutions  of  some 
of  the  problems  that  have  evolved  out  of  the  general  theme 
of  the  measure  of  a  random  set* 

2*1*  Stevens*  Solution  of  the  Problem  of  Coverage  of  a 
Circle  of  Unit  Circumference  • 

W*L*  Stevens  (1939)  considered  n+1  points  chosen 
independently  and  at  random  on  the  circle  of  unit  circumference. 
To  each  point  an  arc  of  length  D  (0  <  D  <  1)  was  attached, 
the  point  being  at  the  clockwise  end  of  the  arc  so  that 
distances  were  measured  anti-clockwise  around  the  circle*  He 
then  posed  the  question:  what  is  the  probability  that  every 
point  of  the  circle  is  included  in  at  least  one  of  the  arcs , 
i*e*  what  is  the  probability  that  there  does  not  occur  a  gap 
after  any  of  the  arcs  (or,  in  terms  of  random  sets,  what  is 
the  probability  that  the  random  linear  set  formed  by  taking 
the  intersection  of  the  point  set  sum  of  the  n+1  arcs  with 
the  unit  circle  is  equal  to  one)  ?  He  first  proved 


-  11  - 


^ ♦  Suppose  that  n+1  events  (which  are  not  independ¬ 
ent)  can  each  occur  in  two  ways,  denoted  by  H  and  T  .  The 
probability  of  exactly  h  H*s  and  n-h+1  T»s  is 

(2.1.1)  (nJ;1){f(h)-(n+J-h)  f(h+l)+  (n+l-h)  f(h+2)+  ...+(-l)n+1f"(n+l)} 

where  the  probability  that  any  h  specified  events  are  H  , 
whatever  be  the  remainder,  is  known  to  be  f(h)  (h=  0,1, . • *,n+l) . 

(Note  that  (2*1«1)  may  be  written  as 

%n+l-h 

h  P  ' 

where,  after  expansion,  we  make  the  substitution  pr=  f (r)  )  • 

Proof :  Let  f(h,t)  (h,t  = 0,1, • * • ,n-f-l  \  h+ t  <  n+1)  denote 

the  probability  that  any  h  specified  events  are  H  ,  any 
other  t  specified  events  are  T  ,  and  the  remaining  are 
either  H  or  T  •  First,  Stevens  evaluated  the  probability 
f(h,l)  of  h  specified  H?s  and  one  specified  T  ♦  The 
admissible  configurations  are  those  with  the  h  specified 
H’s  ,  excluding  those  which  have  an  H  at  the  event  which  is 
to  be  a  T  ,  i^e*  excluding  configurations  with  h+1 
specified  Hfs  ♦ 


Hence 


•  :  '  l  ‘  -  H  .  -1  '■  • 


Iff  1  + 


' 


- 


- 
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f(h,l)  =  f(h)  -  f(h+l) 
-  -  A  f (h)  . 


Using  a  similar  argument. 


f(h,2)  =  f(h,l)  -  f (h+1,1) 
=  -4f(h,l) 


=  A2f(h)  . 


Generalizing,  he  obtained 

f(h,t)  =  (-A)*  f(h) 

=  (l-E)*  f(h) 

(2.1.2)  =  f(h)  -(  *)  f(h+l)  +  (*)  f(h+2)  -  ...  +  (-l)tf(bft)  . 

Assuming  (2*1*2)  to  be  true  for  all  h  and  a  particular  t  , 
it  can  be  proved  easily  by  induction  that  (2*1*2)  is  also 
true  for  all  t  •  Hence,  the  probability  that  h  specified 
events  are  H  ,  and  the  remainder  are  all  T  ,  is 

n+l-h  n+l-h 

(2.1.3)  f(h, n+l-h)  =  f(h)  -  (  1  )f(h+l)+  (  2  )f(h+2)-... 


+  (-1) 


n+l-h 


f(n+l) 


* 
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Now,  the  probability  of  exactly  h  unspecified  H»s  and  the 

remainder  all  T  (and  hence  (2.1,1))  is  obtained  directly  from 

equation  (2,1.3)  by  noting  that  there  are  (  )  ways  of 

h 

selecting  h  events  from  n+1  distinguishable  events,  which 
completes  the  proof. 

Next ,  Stevens  considered  the  existence  or  non¬ 
existence  of  a  gap  after  any  specified  arc  to  be  the  alterna¬ 
tives  H  and  T  respectively .  On  making  this  identification, 
it  follows  that  the  probability  that  a  gap  does  not  occur 
after  any  arc  is  given  by  (putting  h=0  into  (2,1.1)) 

(2.1.4)  f(0)  -(n^)f(l)  +  (n+3)f(2)  -  ...  +(-l)n+1f(n+l)  , 

where  f(i)  is  the  probability  that  there  are  gaps  in  i 
specified  places  (and  in  the  remaining  places  there  may  or 
may  not  be  a  gap).  Therefore,  in  order  to  obtain  the  required 
probability,  we  need  only  determine  f(i)  ,  Stevens  began  by 
finding  the  probability  f(l)  that  there  is  a  gap  after  the 
r  th  arc.  The  restriction  that  there  shall  occur  a  gap  after 
the  r  th  arc  is  equivalent  to  the  restriction  that  no  arc 
shall  occur  in  the  portion  of  the  circumference  of  length  D 
terminating  at  1  .  For  any  configuration  obeying  the  former 
restriction  may  be  converted  to  one  obeying  the  latter 
restriction  by  moving  arcs  r+l,r+2, . . . ,  n+1  clockwise  through 


' 
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a  distance  D  „  Conversely,  the  latter  may  be  converted  to 
the  former  by  moving  arcs  r+l,r+2, « * c ,n+l  anti-clockwise 
through  a  distance  D  „  The  required  probability  is  there¬ 
fore  the  probability  that  the  n  arcs,  excluding  the  first, 
shall  fall  in  a  portion  of  the  circumference  of  length  1-D  , 
i«e0 


f(l)  =  (1-D)n  . 

Notice  that,  in  keeping  with  the  definition  of  f(h)  used  in 
the  theorem,  the  possible  configurations  contributing  to  this 
probability  include  those  which  contain  gaps  elsewhere  than 
in  the  specified  place <,  In  the  same  way,  he  found  the 
simultaneous  probability  f(2)  of  gaps  after  r^  and  r 2  , 
where  r2  ^  rl  9  as 

f ( 2 )  =  (1  -  2D)n  . 

(Since  the  restriction  that  there  be  gaps  after  r^  and  r 
is  equivalent  to  the  restriction  that  no  arc  shall  occur  in 
the  portion  of  the  circumference  of  length  2D  terminating 
at  1  ;  for  a  configuration  of  the  former  kind  can  be  converted 
to  one  of  the  latter  by  moving  arcs  ^+1,^  +  2  ,  , .  0  ,r2 
clockwise  through  a  distance  of  D  ,  and  arcs  r2+l,r2+2, . . . ,n+l 
clockwise  through  a  distance  of  2D  ,  or  conversely.  As 
before,  configurations  having  gaps  in  the  specified  places 


*- 


■  ■ 


. 
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and  elsewhere  are  included  in  the  above  probability) •  The 
generalization  is  now  obvious,,  If  there  are  gaps  in  i 
specified  places,  the  distances  between  each  pair  of  arcs 
enclosing  a  gap  may  be  reduced  by  D  ,  to  produce  a 
configuration  with  the  portion  of  length  i D  terminating  at 
1  containing  no  arcs;  and  conversely .  Hence 

f  (1  -  iD)n  i  <  k  , 

(2.1.5)  f(i)  -  j 

i,  0  i  >  k  , 


where  k  is  the  greatest  integer  less  than  l/D  „ 

Therefore,  from  equations  (2.1.4)  and  (2.1.5), 
the  probability  that  every  point  on  the  circumference  of  the 
circle  is  included  in  at  least  one  of  the  arcs  is 


(2.1.6)  1  -  (  n+1) (l~D)n  +  (  n+1) (l~2D)n 

1  2 


+  (-l)k(  n+1)(l-kD)n 
k 


Notice  that  from  equations  (2.1.1)  and  (2.1.5), 
we  may  write  down  immediately  the  frequency  distribution  of 
the  number  of  gaps.  For  example,  the  probability  that  there 
are  exactly  j  gaps  (j  = 0,1, . . . ,k)  is 


n+l-j 


n+l-j 


(2.1.7)  (  )-!  (l-jD)n-(  )(l-(j+l)D)n+  (  )(l-(j  +  2)D)n  +  .. 

j  1  1  2 

+  (-l)k~d(n+1  3)(l-kD)nj 


;  .  - 
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2  o  2  o  Moments  of  the  Measure  a£  .a  Random  Set  „ 

In  the  first  of  two  papers,  H„E,  Robbins  (19440 
derived  formulae  for  the  moments  of  the  measure  of  any  random 
set,  and  applied  these  formulae  to  find  the  mean  and  variance 
of  a  random  point  set  sum  of  intervals  on  the  real  line,, 

These  formulae  were  used  by  Robbins  (1945)  and  others  to 
obtain  more  generalized  results „ 

Robbins9  most  important  contribution  was  the  proof 
of 

Theorem  2„  Let  X  be  a  random  Lebesgue  measurable  subset  of 
Euclidean  n-dimensional  space  En  ,  with  measure  p(X)  .  For 
any  point  x  of  En  let  p(x)  =  Pr(x  £  X)  0  Suppose, 
furthermore  that  the  function 

f  1  for  x  €  X  , 

(2,2,1)  g(x,X)  ={ 

l  0  for  x  f  X 

is  a  measurable  function  of  the  pair  (x,X)  0  Then  the 
expected  value  of  the  measure  of  X  ,  p(X)  ,  will  be  given  by 
the  Lebesgue  integral  of  the  function  p(x)  over  En  ,  i,e, 

E(yi(X))  =  J  p(x)  du(x)  . 

En 


(2.2,2) 


Iv 


*  )yb  1 
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Proof :  Denote  by  T  the  space  of  all  possible  values  of  X 

and  let  there  be  defined  a  probability  measure  p(X)  on  T 
so  that,  for  every  p-measurable  subset  S  of  T  ,  the 
probability  that  X  shall  belong  to  S  is  given  by  the 
Lebesgue  -  Stieltjes  integral 

(2.2.3)  Pr(XtS)  =  J  C  (X)dp(X)  , 

T  * 


where  the  integrand  is  the  characteristic  function  of  S  , 


(2.2.4)  C  (X)  =- 

o 

V 

For  every  point  x 
as  the  set  of  all  X  in  T 
Then  for  every  point  x  in 

(2.2.4) 


1 

for 

X  6  S  , 

0 

for 

X  {  S  . 

of 

E  , 
n  9 

Robbins  defined  S(x) 

which  contain  the  point  x  . 

E  we  have  from  (  2  *  2  „  l)  and 
n 


(2.2.5)  g(x,X)  =  Cg(x)(X)  . 


Denoting  the  Lebesgue  measure  in  En  of  the  set  X  by 
p(x)  ,  we  have 


. 
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(Zo2  06)  E(n(X))  =  J*  u(x)dp(X), 

T 

and,  from  (2„2.1),  we  have 

( 2  o  2  o  7)  |i(X)  =  J  d|j(x)  =  J  g(x,X)dll(x)  . 

X  E 

n 

Now,  assuming  that  the  function  g(x,X)  is  a  X  p  - 

measurable  function  of  the  pair  (x,X)  in  the  product  space 

of  E  with  T  ,  it  follows  from  Fubinivs  theorem  and  from 
n 

(2o2  06)  and  (2o2<,7)  that 

E(p(X))  =  J  /  g(x,X)dp(x)dp(X) 

T  En 


(  2  o 2  0  8  )  =  J*  J  g(x,X)dp(X)dp(x)  =  J  g(x,X)dtip(x,X)  # 

E  T  EXT 

n  n 

But,  from  (2o20l)  and  (20203), 


( 2  o 2  o 9 )  J  g(x,X)dp(X)  = Pr(X 6  S(x) )  =  Pr(x €  X)  , 

T 

and  hence  (2<,2<>2)  follows,  ic.e0 


. 
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E(p(X))  —  J*  Pr (x € X)dp(x)  • 

E 

n 

Robbins  next  generalized  the  above  fundamental 
result  to  obtain  similar  expressions  for  the  higher  moments 
of  u(X)  .  The  formula  for  the  m  th  moment  of  p(X)  was 
given  as 


(2,2.10) 


E(p  (X))  J  p(x1>x2, o . . ,xm)  dp  (x1,x2 . . . ,xm)  ^ 
E 

mn 


where  ,x0  ,  . . .  ,xj  denotes  Lebesgue  measure  in  E 

±  2 '  *  m  “  mn 

and  where 


p(x± 


900 


x  )  =  Pr (xn  6  X  and  x0  6  X  ...  and  x  E  X) 


m 


m 


We  are  now  in  a  position  to  apply  the  above  results 
of  Robbins  to  the  random  set  X  described  at  the  beginning 
of  Chapter  1  .  Recall  that  we  can  define  this  random  set  X 
as  follows:  n+1  points  are  chosen  independently  and  at 
random  on  the  circle  of  unit  circumference.  To  each  point  an 
arc  of  length  D  (  0  <  D  <  l)  is  attached,  the  point  being  at 
the  clockwise  end  of  the  arc.  The  set  X  is  that  portion  of 
the  circle  of  unit  circumference  that  is  covered  by  the 


■ 
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n+1  arcs.  Now,  if  a  set  is  measurable,  then  its  characteris¬ 
tic  function  is  a  measurable  function 0  Hence,  since  the  set 
X  defined  above  is  a  measurable  set,  the  conditions  of 
Theorem  2  are  clearly  satisfied,,  The  probability  that  any 
point  x  on  the  circle  of  unit  circumference  shall  be 
contained  in  the  i  th  arc  is  clearly  D  0  Hence 

(2„2oll)  Pr(x €  X)  =  Pr(x  belongs  to  at  least  one  of  the  n+1  arcs) 


=  1  -  Pr(x  does  not  belong  to  any  of  the  n+1  arcs) 


=  1-  (l~D)n+1 


9 


so  that,  for  any  point  x  on  the  circle  of  unit  circumference. 


Pr(x  €  X)  =  1-  (1-D)n+1 


( 2  o  2 „12 ) 


o 


From  (2,2,2)  it  follows  that 


( 2  o  2 „ 13) 


0 


To  evaluate  E(p  (X))  we  make  use  of  the  identity 


( 2  o 2  o 14)  Pr ( A  and  B)  =  Pr(A)  +  Pr(B)  +  Pr(neither  A  nor  B)  -  1 


9 


c 
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which  holds  for  any  two  events  A  and  B  0  It  follows  from 
(202ol2)  and  (2020l4)  that  if  x  and  y  are  any  two  points 
on  the  circle  of  unit  circumference,  then 

(2o20l5)  p  ( x ,  y )  =  Pr(x€X  and  y  €  X ) 

=  Pr(x£X)  +  Pr(y€X)  +  Pr(xfX  and  y  f[  X)  -  1 

=  1-  2(l~D)n+1+  Pr (x  <f  I  and  j  j  X)  . 


Let 

( 2  o  2  o  16 )  h(x,y)  =  Pr(x|  X  and  y  X)  , 

and  assume,  without  loss  of  generality,  that  y  >  x  <>  From 
(2o2ol0),  (2„2ol5),  and  (2„2ol6),  we  want  to  evaluate 


1  1 

( 2  o  2  o  1?)  E(vjl2(X))  =  J*  J  p(x,y )  dy  dx 

0  0 

=  J  J  [l-2(l-D)n+1  +  h(x,y)]  dy  dx 

0  0 


=  1-2(1-D)n+1+  2 


1  1 

J  J  h(x,y)  dy  dx 
0  x 


o 


For  D  <  y 


it  may  be  shown  that 


. 
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( 2  o  2  o  1 8 ) 


(l-(D  +  y-  x)  ) 


n+1 


=  J 


(1-2D) 


n+1 


V. 


( 1-  (  D  +  x  +l-y  )  ) 


n+1 


for 

y  -x 

<  D 

for 

D  <  y  -  x 

<  1  -  D 

for 

y  -  x 

>  1  -  D 

Hence 9  from  (202017),  when  the  latter  integral  in  that 
expression  is  evaluated,  we  obtain 


E(u2(X))  =  1-2(1-D)n+1+  (  2(l-D)n+2  +  n(l-2D)n+2)  (D  <1)  . 

Combining  this  with  (202o13),  we  obtain  the  variance  of  p(x) 
as 


(2.2.19)  Var  (W(X))  =  E(u2(X))  -  E2(y(X)) 


1 

n+2 


(2(1»D)I1+2  +  n(l-2B)n+2 


(1-D)2n+2  (D  < 


If  D  >  ^ 


it  can  be  shown  that 


(l-(D+y-x 

))n+1 

for 

y  -  x  <  1 

( 2  o  2  0  20 )  h(x,y ) 

0 

for  1  - 

-  D  <  y  -  x  <  D 

L  ( 1  ~  (  B  +  x  +1 

-y))n+1 

for 

y  -  x  *>  D 

from  which,  together  with  (2o20l/)$  we  obtain 


*  ri 


((X)^q)3  -  ((X)q)  n*V 


23  - 


E(u2(X))  =  1-2(1-D)n+1+  2(1-D)n+2  (D  >  A  )  , 


and 


(2.2.21)  Var  (p(X))  =  ^  (1-D)n+2  -  (i_D)2n+2  (D  >  1  )  # 


2.3  "The  Probability  Distribution  of  the  Measure  of  &  Random 
Linear  Set" 

Consider  now  the  problem  of  determining  the 
probability  function  of  the  measure  of  a  one-dimensional 
random  set.  Let  n  points  z^(i=l,2>...,n)  be  selected 
independently  and  at  random  from  the  interval  (0,1)  ,  the 
distribution  of  any  being  the  uniform  distribution  with 

distribution  function  given  by  (l.l),  i.e. 

F(z)  =  z  (0  <  z  <  1) 

Let  these  n  points  be  arranged  in  increasing  order  of 
magnitude  as  xi(i=  1,  2,  » .  . ,n)  ,  and  with  each  x±  let  there 
be  associated  an  interval  of  length  D(0  ^  D  ^  1 )  centered 
at  xi  .  Let  p(X)  denote  the  measure  of  the  set  X 
consisting  of  the  point  set  sum  of  the  intervals 

( x .  —  7T  D  .  x.+  tD)  (0  <  D  <  1  ;  i-1,2,...  ,n)  • 

l  2  l  2 


« 


We  now  present  in  detail  a  paper  by  D^F,  Votaw  Jr„  (1946),  in 
which  he  calculated  the  probability  function  of  the  random 
variable  p(X)  „  This  is  a  very  useful  and  important  paper 
in  that  the  general  method  used  by  Votaw  in  this  paper  leads 
us  to  the  solution  of  the  similar  one-dimensional  but 
circular  problem  that  was  described  in  the  introductory 
chapter  0 

Votaw  began  by  making  the  transformation 

=  xi 

yi  =  xi+l_xi  (i- 1,2,..., n-1)  , 

and  letting  y=p(X)-D  0  The  range  of  y  becomes  0  <y  <m  , 
where  m  denotes  the  minimum  of  1  and  (n-l)D  „  It  is 
clear  that  y  may  be  expressed  as 

n-1 

( 2  03  el)  y  =  I  m(y , ,D)  , 

i=l 

where  m(yisD)  denotes  the  minimum  of  y±  and  D  0  From 
the  theory  of  order  statistics  it  can  be  shown  that  the  joint 
density  function  of  x-^ ,  x^  9  o  <>  o  9  xr  is  nlCo^x^^x^^ooo^  xn<*  ^  ^ 

and  since  3(yoSy^>  «  °  °  >yn_l^/^^xlsX2s' °  °  oXn^  ~  1  *  it  is  clear 

i  is  n! 
n-1 


that  the  joint  density  function  of  yo,y1Sooo,y 


. 


,  v)  to  nvfi* 
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If  m  -  (n~l)D  ,  it  is  clear  that  there  will  be  a 


finite  probability  that  the  set  X  will  consist  of  the  union 
of  n  disjoint  intervals  of  length  D  and  hence,  the 
probability  function  of  y  will  have  a  saltus  point  at 
y=(n-l)D  .  Now,  given  that  m  =  (n-1 ) D  ,  y  =  (n-l)D  if  and 
only  if  y^  >  D  (i = 1,2, . . . ,n-l)  .  In  order  to  determine  the 
amount  of  the  saltus  at  y*  (n-l)D  ,  Votaw  first  fixed  yQ 
and  calculated  the  volume  of  the  (n-l)  -  dimensional  region 
in  which  any  point  (y  ,  «.  o  o  *yn_i_)  satisfies  the  condition 

y^  >  D  ,  y^  >  D,00o,  yn  ^  ^  D  .  Recalling  that  we  must  have 

n-l 

E  y.  <  1  ,  this  region  is  clearly  R  in  the  figure  below, 
i=0  1 

and  its  volume,  V  say  ,  is  given  by 


“  J  J  -oo  J  dyn-1 ° o  ody2dy1 

0  0  0 


plane  y..+y0+y,  =  1  -y 


(i-yo>°>°) 


1*2*3 


*1 


o 


is  tfnioq  euilAt  a  lli  f 
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y  .  —  D 

Now,  making  the  substitution  z.  =  Z - - — r —  (i=l  2  n-l'i 

1  l-yo-(n-l)D  V1  x9*9°ooSn  i) 


n~l 


so  that 


£  z.  <1  (z  >  0,  i  =  l,2 
i=l  1  1 


£  o  o  o  y 


n-l)  ,  we  have 


1  1-z. 


v  =  (l-yo-(n-l)D)n'1  J*  J* 

0  0 


l-Zl-Z2-o..-Zn_2 

of  °  °  >  „  dz , 

0  L  1 


which,  by  Dirichlets  integral,  becomes 


V  = 


(l-yo  -  (n  -1)  D)n_1 
(n-1) ’ 


By  letting  yQ  vary,  it  is  clear  from  the  figure  above  that 


1-  (n~l)D 

(2.3.2)  Pr(y=(n-1)D)  =  J  ( l-yQ  -  (n-1 )  D)n-1dyo 

0 

=  (l-(n-l)D)n  (n-l)D  <  1  . 


Now,  the  sample  space  on  which  the  random  variable 
y  is  defined  consists  of  all  points  (yQ syl 9  0  °  ° ,yn-l ^ 
belonging  to  the  n-dimensional  tetrahedron  defined  by 

n~l 

y±  >  0  (i=  0, 1 , 0  .  . ,n-l)  and  X  y.  <  1  0  As  we  shall  see 

i— 0  1 


?  ■  *  <  v: 
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later,  the  probability  that  Y<y<Y+AY  (where  Y  <  m 
and  AY  denotes  an  arbitrary  small  positive  increment  in  Y) 
can  be  evaluated  by  determining  volumes  of  certain  regions 
contained  in  this  tetrahedron.  To  provide  a  starting  point 
for  the  derivation  of  the  probability  function  of  y  , 

Votaw  considered  the  following  conditions. 

(a)  qD<Y<(q  +  l)D  (q=0,l,...,M;M  denotes  the 

minimum  of  (n-2)  and  the  greatest 
integer  less  than  l/D  )  . 


For  any  particular  value  of  q  ,  he  decomposed  the  set 
y  ,Yl, . . « >yn-i  into  three  subsets  (yQ)  ,  (y± ,y2 , . . . ,Y j , 

and  (yj+1  yd+2J'oosyn-l)  such  that  with  yo  fixed 

(k)  yu^D  (u  =  l,2,,.,sj  ;  j  <  q)  , 

(c)  yy  <  D  ( v = j+l,o+2, . o o ,n-l) 

n-1 

Now,  since  E  y .  <  1  ,  we  have 

i=0 


3 

E  y 

u=0 


u 


+ 


n-1 

E  y  <  1 

V=j+1  V 


9 


' 
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and,  as  a  consequence  of  (b)  and  (c)  we  have,  from 
equation  (2o30l) 


n-1 

y  =  3 D  +  E  y  ; 
v=j+ 1 


thus 


3 

z  yu+y  -Jd  <  i  * 

u=0 

From  this,  we  have 

3 

( d  )  Xy  <l-y-y  +  jD  0 

u— 1  ° 

Now,  let  A. (y,y  )  denote  the  volume  of  the  j -dimensional 
region  in  which  any  point  (y-L,y2  > <> <>  o  ,y^  )  satisfies 
conditions  (b)  and  (d)  ,  and  let  B  9(y)  denote  the 

volume  of  the  (n-j -2) -  dimen sional  region  which  is  the 

n-1 

intersection  of  the  hyperplane  X  y  =  y- j  D  with  the 

v=j+l  V 

(n-j-l)  -  dimensional  cube  given  by  condition  (c)  0 

Let  us  first  evaluate  A.  (y,y  )  <>  For  convenience, 

vl  V 

let  h  =  l-yQ~y +  j D  °  Then  from  conditions  (b)  and  (d)  , 
this  volume  is  simply  (see  derivation  of  (20302)) 


I  ♦  £ ) 


•  . 


. 
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h  h-yi  h-y1~y2-  *  .  „  -  y_. 

J*  J*  *  °  •  J  dy^. .  „  ,dy2dy1 


Letting  y^  h§^  >y2  ~  ^2  9  0  9  9  >yj  ~  h?j  »  we  have 

1  1_  ^1  1_^L~  ^2~  °  °  ° 

Aj(y»y0)  =h1  J*  J...  J  d?  d?2d?1  . 

D/h  D/h  D/h 


•  1 

Now,  on  making  the  substitution  z.  =  — 1 — —  We  have 

__  * 

1  “J  h 

o 

S  z .  <  1  and  z .  >  0  (1*1,2,. .,,3),  so  that 
i=l  1  1 


1  1-z-j^  1 


Ad(y,y0)  =  hJ(1-  J  J  ... 

00  0 


l’z2~°”~z3-l 


d z ^  o  o  ©  d  Z‘2 dz-j 


(2  03o3) 


(h-jD)^1 


(1-y0-y)3 


a! 


Note  that  the  range  of  yQ  is  0  <  yQ  <  1-y  •  To  evaluate 

B  _._9(y)  ,  Votaw  used  a  result  from  a  paper  by  P*  Hall  (1927), 
n~j  z 

in  which  it  is  proved  that  the  volume  cut  from  the  hyper plane 


■:  r;i 


30 


N 

T.  x  =p  by  the  unit  cube  0  <  x.  <1  (i  =  1, 2 , „  . .  ,N )  is 
i=l  1  l 

given  by 


(2.3.4) 


v(p) 


(N-i): 


E  (-l)r  (  "  )  (p-r)N_1  , 

r=0  r 


where  K  is  equal  to  the  greatest  integer  less  than  p  . 
Now,  we  wish  to  determine  the  volume  cut  from  the  hyper plane 


n-1 

s  yv  =  y-dD 

v=j+l  v 


by  the  hypercube  0  <  yv 
y  =z  .  ( v=j+l, « . o ,n-l) 

V  ’  J 

the  hyper plane 


<D  (v=j  +  l,  .  ..,n-l)  o  Let 

* 

,  then  we  want  the  volume  cut  from 


n-j-1 

E  z .  =  y-jD 
i=l  1 


by  the  hypercube  0  —  ®  (i~l,2, 0  0  <»  ,n-j  l)  <,  In  order  to 

apply  Hall’s  result,  let  us  transform  the  above  hypercube  to 
the  unit  hypercube  by  letting  =  z^/D  *  Since  0  <  D  <  1 
this  transformation  is  an  expansion,  so  we  can  find  the 
volume  cut  from  the  (n-j -2) -dimensional  hyper plane 


y 


U.C.s) 


» 
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n-j-1 


E  x.  =  x 


i=l 


l  D 


-  3 


by  the  unit  hypercube  0  <  xi  <1  (i=l,2,  n-j-1 )  using 

J  ey 

Hallfs  result  and  then  multiply  the  result  by  Dn  J  „ 
Hence,  using  Hall,  we  have  p  =  ~  -  j  ,  N=  n-j-1  ,  and  since 


qD  <  y  <  (q+l)D 
q- 3  <  q  "  3  <  ( q+1 )  ”  3 
q-j  <  P  <  (q+1)  -  3  5 


we  have  K  =  q-j  ,  so  that,  from  (2„3<,4) 


- —  q -3  n-j-1  n-j-2 

b  -i_o(y) =  Pn~3~2  'Jn~3—  s(-Dr(  „  )(§-J-r) 


n-j-2 


(n-j-2)I  r=0 


(2.3.5) 


I -  q  3  n-j-l  r  -,n- j-  2 

=  *±i  s(-Dr(  r 

(n-j-2);  r=0 


)  j^y-D(  j+r)  J 


A . (y ,y  )  and  B  .  0(y)  define  contents  of  regions  in 
3  w,,7o  n-j-2 

orthogonal  subspaces  of  the  (n-l) -dimensional  tetrahedron 


n-l 

£  y  .  1  —  y  »  The 

i=l  1 


content  of  the  j  +  (n-j-2)  — n  2 


■ 
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dimensional  region  satisfying  conditions  (b)  ,  (c)  ,  and  (d) 
is  therefore  Aj (y >yQ)  Bn_ j _2(y )  •  Therefore,  the  probability 

that  Y  <  y  <  Y  +  AY  and  that  conditions  (b)  ,  (c)  ,  and 
(d)  are  satisfied  is 

Y+AY  i-y 

(2.3  =  6)  n *  J  Bn_j_2(y)  J  Aj(y,yo)  dyQ  dy  # 

y=Y  yo  =  °  Jn-j-l 


To  illustrate  why  the  factor  1/Vn- j-1  occurs  in  (2,3,6)  , 
consider  the  following  example.  Suppose  we  let 

n  =  7  t  D  =  ^  ,  and  suPPose  we  take  a  particular  case  where 

q=4  and  j  =  3  ,  Hence,  conditions  (a)  through  (d)  become 


(*»)  yl  ^  6  >y2  -  6  >y3  ~  6  , 

(c)  y4  <  6  ’  y5  <  6  '  y6  <  6  > 

(d)  yx  +  y2  +  y3  -  1  ” yo  "  y  +  2  * 


and  we  also  have  y^  +  y^  +  y^  y  2  ®  Now,  B^(y)  is 

equal  to  the  content  of  the  2-dimensional  region  containing 
all  points  (y^^^)  which  satisfy  condition  (c)  and 


- 


<i  9  .)  ->  >C  6  0:4  ct  ©ivoo  q  u  ..  , 
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the  equation  Z  y  y  -  «  >  i*e*  B?(y)  is  the  content  of 

v=4  z 

the  plane  Z  y  y  -  cut  off  by  a  cube  (at  the  origin)  of 

v=4  z 


side 


We  want  to  evaluate  the  probability  that  Y  <  y  <  Y  +  AY  and 
that  conditions  (a)  through  (d)  are  satisfied.  This  is 
given  by  equation  (2.3.6)  which  involves  the  volume  between 

6  i  6  x 

the  planes  Z  y  =y-  “  and  Z  y  =  y+dy  -  ”  .  To  determine 

I  V  Z  _ t  V  z 

v=4  v=4 

this  volume,  we  need  only  multiply  B2(y)  by  the  perpendicular 
distance  between  the  planes  (this  is  assuming,  of  course,  that 
dy  is  small  enough  that  we  can  neglect  the  volume  gained  by 


. 


y  fit  :q  ©i  ‘J 


ei>  I”  tfo  nso  sw  rict  rgi  ri  t  bma  ai  y 
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not  considering  the  effect  of  the  boundary  of  this  region  with 
the  cube).  This  distance,  dz  say,  is  simply 


Hence,  the  volume  of  the  region  is 


B2(y)dz—  B^Cy)  ^  .  In  general  l/Vn-j-l  is  the  cosine 

of  the  angle  between  a  coordinate  axis  and  the  normal  to  the 
hyperplane  y  .+^ +  . . . + ynl  =  y-jD  •  If  any  other  of  the 

(  n  ^  )  combinations  of  the  j  yfs  had  been  specified  in 
3 

(b)  and  the  (n-j-l)  complementary  y*s  had  been  specified 
in  (c)  ,  the  corresponding  A_.(y,yQ)  and  Bn_j_2(y)  would 

be  equal  to  those  in  (2.3.6)  ,  so  that 


Pr(Y  <  y  <  Y  +  AY) =  n» 


q 

X  ( 

3=0 


Y+  Ay  l-y 

nj1}  J  Bn-j-2(y)  J'  Aj(y’yo)dyo 

y=Y  y  =  0 

J  o 


for  qD  <  Y  <  ( q+1  )  D  ,  Y  ^  m  ,  ( q  —  0,1,...,  M)  . 


Therefore,  from  (2.3.3)  and  (2.3.5)  we  have 


Y  +  Ay 


Pr(Y  <  y  <  Y  + 


q  .  „  1  q-3  .  _  r  -|n”3"‘2 

AY)=n;  Z  (n.1)J  (-.-i-TTt  S  (-Dr(n  J-1ly-D(d+rM 

j=0  3  Y  1  '*  2J'r=0  r  J 


i-y  -j 

(l-y  -y) 


.  / 


0 


P  "  dy_  dy 
i  t  O 


Ci 


■  •  ,  '  (  ) 
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q  q-o’ 

n  £  I  (-i)r  (  n"l)(  n-d-l  x/  n-1 


j=0  r=o 


)(  tt_A) 

r  M  o  +1 


Y+*Y  n  ^  7 

.  J  [y-D(d+r)]n  d  Vy^dy  . 
Y 


Applying  the  mean  value  theorem,  there  exists  a  ?  such  that 


Pr(Y<y<Xt^_l=n  l  ,id(.1)r(n-l)(n-d-l)(„;l) 

iY  j=l  r=o  J  r  J+1 


r  -.n-j-2  j+1 

.  [5-D(d+r>]  (1-1) 


where  Y  <  ?  <  Y+  AY  #  Letting  A  Y  -  0  ,  we  have  the 
probability  function  of  y  as 


q  q-j 

(2.3.7)  f  (y) =  n  S  r(-l)r(  )( 

j=l  r=0 


n 


r  (  “-lw  n- j  -1  ^  ^  n-1  ^ 


3+1 


r  -»n-o-2 

j_y-D(j+r)J  (l-y)J+  , 


qD  <  y  <  (q+l)D  ,  (q = 0 ,1, . • . ,M)  ,  y<m. 


Votaw  next  noted  that  if  (n-l) D  <  1  ,  then  fn(y) 
is  not  defined  at  y=(n-l)D  (hence  the  range  of  definition 
of  f  (y)  is  given  as  0  <  y  <  m  ) .  This  is  due  to  the  fact 


. 
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that  there  is  a  saltus  of  amount  (l-(n-l)D)n  at  y=  (n-l)D 
(see  (2<,3o2))o  If  m  =  l  ,  the  range  of  y  in  (2„3o7)  is 
0  <  y  <  1  o  He  then  stated  that  the  probability  function 
f  (y)  is  continuous  over  the  range  0  <  y  <  m  with  the 
exception  of  a  simple  discontinuity  at  y=  (n-2)D  in  the 
case  n  >  3  and  (n-2)D  <  1  .  The  distribution  function  of 
y  is  continuous  with  the  exception,  in  the  case  (n-l)D  <1  , 
of  a  discontinuity  corresponding  to  the  saltus  of  amount 
(l-(n-l) D)n  at  y  = (n-1 ) D 

The  expected  value  of  y  (and  hence  of  p(X))  can 
be  derived  in  the  usual  way  from  f  (y)  or  by  the  use  a 
result  of  Robbins o  Here,  we  present  the  latter  method,. 
Recall  that  the  expected  value  of  the  measure  of  a  random 
Lebesgue  measurable  subset  X  of  Euclidean  n-dimensional 
space  ,  is  given  by  the  Lebesgue  integral 

E(p(X) )  =  J  p(x)  dx  , 

En 

where,  for  any  point  x  of  E  ,  p(x) =  Pr(x  6  X)  ,  provided 


the  function 


g( x,X) 


r  1  for  x  €  X 
w  0  for  x  $  X 


$ 


* 


<  " ; 


"Of 


-  '  ■ 
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is  a  measurable  function  of  the  pair  (x,X)  0  Let  x  be 
any  point  belonging  to  the  interval  (-D/2  ,  1+  D/2)  ,  and 
let  p  ,  say,  be  the  probability  that  x  shall  be  contained 


in  the  i  th 

subinterval  of 

length  D  centered  at 

x^ (i  1,2,oq0, 

n-l)o  It  is  clear,  then,  that 

(  x  +  f 

x  €  (-  2  ^  0  ) 

x  +  2 

x  €  (  0  ,  ;~  ) 

’O 

II 

D 

x  €  (  2  ^  1  “  f  ^ 

1  “X  +  f 

x  €  (  1-“,1  ) 

i  -  +  £ 

l  1  X  +  2 

x  t  (  1  9  1  +  2  )  9 

or 

r  x  +  2 

-  £  <x 

2  —  2 

p  =  i 

D 

£  <r  £ 

2  —  X  <1  "2 

l 

l  1-  x  +  § 

1  -  f  <x  <1  +  f  t 

Now,  from  (2. 

2  o ll)  ,  p(x)  = 

=  Pr (x  €  x)  =  1  -  (l-p)n  ,  hence 

for  0  <  D  <  1 

D/2 

[l-(l-x-D/2)n] 

l-D/2  l+D/2 

E(p(X))  =  J 

dx+J  ^1-(1-D)n  J  dx  +  /  |_1- 

-D/2 

D/2  l-D/2 

D/2)" 


. 
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( 2  o  3 , 8 ) 


=  D  + 


n-1 

n+1 


(1-(1-D)n+1  ) 


o 


2jt  Random  Intervals  on  n  Line 

C0  Domb  (1947)  considered  a  sequence  of  random 
events  occuring  as  points  on  the  line  t  from  -  00  to  +  00  8 
On  assuming  that  the  probability  of  the  occurrence  of  two  or 
more  of  these  points  in  a  very  short  interval  (t  ,  t  +  dt ) 
is  X  dt  ,  one  can  show  that  the  probability  of  finding 
exactly  j  points  (events)  within  a  fixed  interval  of  length 
P  say,  is  given  by  the  Poisson  distribution,  with  parameter 
X  $  x  q  e  q 


p( 3  ,  X  P  ) 


e~XP(XB)J 


Domb  considered  the  case  where  each  event  (or  occurrence  of 
a  point)  consists  of  an  interval  of  length  D  (a  positive 
constant);  an  event  being  characterized  by  its  first  point  * 

He  let  W(x  ,  P)  be  the  differential  coefficient  (with  respect 
to  x  )  of  the  function  T(x,  P)  defined  to  be  equal  to  the 
probability  that  the  covered  portion  of  the  interval  (0  ,  p) 
is  less  than  or  equal  to  x  .  He  then  obtained  the  function 


G(p,z)  =  e^z  &  (w(x,P)} 


(2.4.1) 


9 


where  z  =  0-x  and  £{f(x,z)}  =  G(p,z)  denotes  the 
p-multiplied  Laplace  transform  of  f(x,z)  ,  i.e 

00 

(2.4.2)  P  |  e  px  f(x,z)  dx  =  G(p,z)  , 

0 

as  the  solution  of  an  integral  equation  of  the  convolution 
type.  From  this  solution  he  was  able  to  derive  the  following 
results.  The  probability  that  the  total  interval  (0 ,  0)  is 
covered  is  given  by  the  function 

(2.4.3)  z(p)  =  l-e'XD(l+XP)+e"2XD  [x(p-D)+  -  ..  . 

,  ltr+l  -(r+l)XD  r  r  (g-rD)r+  r+1  (P-rD)r+1 -i 

1  6  lX  r*  X  (r+l)  t  J+* 

all  terms  beyond  0-nD  being  ignored  if  nD  <  £  <  (n+l)D  , 

The  probability  of  the  interval  being  completely  uncovered  is 

e  ^  ^  while  the  probability  that  r  non-overlapping 

events  occur  in  the  interval  (including  no  overlapping  at  the 
end-points)  is 


(2.4.4) 


Xr(£-rD)r 


-X(0+D) 


(0  >  rD)  . 


■ 
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Next,  he  derived  an  expansion  of  W(x,£)  according 
to  the  number  of  events  in  the  interval  (0  ,  £)  ;  the  complete 
solution  of  which,  he  admitted  as  very  difficult*  However,  he 
succeeded  in  deriving  an  expansion  which  gives  the  probability 
that  given  n+1  events  occur  in  ( ~D ,  p )  ,  the  interval 

(0  ,  p)  is  completely  covered.  Using  this  expression  he  was 
able  to  obtain  a  result  of  Stevens,  namely  (2,1,6)  , 


•  •'  -• '  d*  ■ 


CHAPTER  III 


MAIN  RESULTS 

Consider  n+1  points  zi(i=0,l5O„o,n)  chosen 
independently  and  at  random  from  the  interval  (0,1)  ,  the 
distribtion  of  any  being  the  uniform  distribution  with 

distribution  function  given  by  (l„l)  0  These  points  are 
then  arranged  in  increasing  order  of  magnitude  as 
x±(i  =  0,1,  «,  e  .  ,n)  «,  To  each  point  ,  an  interval  I  is 

associated  where  I  is  given  by  (l«2)  0  Let  X  denote 
the  random  set  which  is  the  point  set  sum  of  the  n+1 
intervals  {l }  3  As  was  mentioned  in  Chapter  I,  this 
representation  is  equivalent  to  the  circular  representation 
of  Stevens,  and,  it  is  this  circular  nature  that  will  allow 
us  to  assume,  without  loss  of  generality,  that  xq = 0 
It  is  also  clear  that  it  is  only  the  cyclical  or  overlapping 
nature  of  this  representation  that  differs  from  the  one¬ 
dimensional  non-circular  problem  considered  by  Votaw  0 

This  chapter  contains  the  main  result  of  this 
thesis;  namely  that  of  determining  (using  the  methods  of  Votaw) 
the  probability  function  of  the  measure,  p(X)  ,  of  the  random 
set  X  *  Before  we  derive  this  function,  however,  we  consider 


a  few  related  problems 


- 


■ 


' 
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3 0.1  Expected  Value  of  the  Measure  q£.  &  Random  Linear  Set 

and  Related  Problems  . 


First,  using  methods  different  from  those  of 
Robbins,  we  again  derive  equation  (2.2*13)  ,  i.e.  we 
calculate  the  expected  value  of  p(X)  .  It  is  clear  that 
p(X)  may  be  expressed  as 

n  -1 

(3.1.1)  u(x)  =  c(Xl)  +  c(xi+1-Xi)  +  c(l-xn)  , 
where 


(3.1.2) 


0  <  x  <  D  , 


D  <  x  <  1  . 


Hence,  from  (3.1.1)  ,  we  have  the  expected  value  of  p(X) 

as 

n-1 

(3.1.3)  E(|jl(X))  =  E(c(x1))  +  ^  E(c(xi+1-  xi))  +  E(c(l-xn))  . 

Now  z-*  z  .  .  .  z  is  a  random  sample  from  a  continuous 

population  (xq  was  fixed  zero)  ,  each  being  uniformly 

distributed  over  the  interval  (0,1)  .  The  r  th  largest  of 

these  sample  values,  namely  xr  ,  is  the  r  th  order 

statistic.  It  can  be  shown  that  the  density  function  gr(xp) 

of  x  is  given  by 
r 


* 
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r-1  oo 


(301o4)  g  (x  ) 
r  r 


fr-l)”fn-r) '  l  I  f(z)dz]  f(zr4  f  f(z)dzJ 

0  '  3  —00  “*  Xf 


n-r 


for  <  00  s  where  f(z)  is  the  common  uniform  density 

function  over  (0,1)  of  the  random  variables  z.  (i  =  1,2,00^)  0 

i 

It  can  also  be  proved  that  the  joint  density  function  of 
for  r  <  s  is  given  by 


x 


(30105)  g(x  ,x  )  = 

r  s 


n«f  (zr)f(zs)  J-  r 

(r-1) ® (s-r-1) *(n-s) J  L  Jrof(z)dz  1 


r-1 


xs  s-r-1  oo 

|"  J  f(z)dzj  [  J  f(z)dzj 


n-s 


for  —oo  <  x  \  x  '  -  c 
r  s 


<  x  <  00  0  It  is  clear,  therefore,  that 


E(c(x-l))  J*  c(x-.  ) g-.  (x-j  ) dR 
R 


where  R  is  the  region  in  which  the  density  function 
g  (x  )  is  defined ,  From  (30lo2)  and  (30lo4)  we  obtain 


D  1 

E(c(x,))  =  n  |  x1(l-x1)“  1dx1  +  nD  J  (l-x^11  dxx 

0  0 


(3=1.6) 


n+1 


(1-D) 

n+1 


n+1 


■ 
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Similarly ,  it  may  be  shown  that 


(3.1.7) 


E(c(l-xn)) 


_L-  _  (1-D)n+1 
n+1  “  n+1 


Now,  consider 


(3.1.8) 


E(c(xi+1~  x±)) 


-  J  J  c(. 


1+1 


xi ) 


g( 


bi+l 


x±)dR/ 


R 


where,  from  (3,1,2) 


(3,1,9) 


c(xi+l  xi  ^ 


'  xi+l  "  xi 


D 


0  <  x 
D  <  x 


i+1 

i+1 


<  D 

<  1  , 


and,  from  (3,1,5) 


(3,1,10)  g(  x 


i+1 


x±) 


i-1 


n 


( i-1 ) l( n-i-1 ) 


(1 


xi+i } 


n-i-l 


0 


<  x .  < 

l 


1  : 


1<X 


and,  where  R^  is  the  triangular  region  of  the  (x. ,x. ,, ) 

i  l+l 

plane  bounded  by  the  lines  x^  =  0  ,  =  1  ,  and  xi+i  =  xi  » 

In  this  region,  however,  the  function  c(x^+-^-x^)  assumes  two 
different  values,  namely  the  constant  value  D  in  the 
triangular  region,  R^  say,  bounded  by  the  lines 


a  >  £  0  1 

' 

. 
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x^-0,x^+^-l  and  x^+^=  x^+ D  ,  and  valu©  xi+l  ”  xi 

in  the  trapezoidal  region,  R*  say,  bounded  by  the  lines 
=  0  ,  x^+^  =  x^  +  D  ,  x^+^  =  D  ,  and  x^+^  =  x±  •  Since  these 
regions  are  disjoint,  we  may  write  (3,1.8)  as 


E(c(xi+l”xi))  =  JJ  c(xi+1“Xi)g(xi+1,xi)dR-[  +  J  j  c(xi+1-xi)g(xi+1,xi)dR2  . 


“i 


R 


Using  (3.1.9)  and  (3.1.10)  ,  this  becomes 


(3.1.11)  E( c(x^+^-x^ ) ) 


1-D  1 


n-i-1 


•  {°  X  X  xi"  (1_xi+i)''  *  *dxi+idxi  + 


+ 


+ 


0  x.  +  D 
1 

D  xi+l 


x  x  ( xi+1-xi )  xi‘ 1  ( i-xi+1  r^dx.  dxi+1+ 


0  0 


1  x 


i+1 


X  X  (xi+i-xi)x-  1(l-xi+1)n'1‘1dx.dx.+1} 


D  x.  ,,-D 
l+l 


Evaluating  the  integrals  of  (3.1.11)  and  combining  with 
(3.1.6)  and  (3.1.7)  we  readily  obtain 


. 
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(3.1.12)  E(u(X))  =  1  -  (1-D)n+1  . 

The  result  given  by  (3.1.12)  suggests  the  follow- 
question •  how  many  points,  n+1  ,  would  be  required  in 
order  to  obtain  an  expected  coverage  of,  say  C  ?  To 
answer  this  question,  we  need  only  determine  an  interger  n 
such  that 


C  =  1  -  (1-D)n+1  . 


It  is  clear  that 

(T.  1  TO  n  =  log  (1-0  ~  log  (1-D) 

n  log  (1-D) 

so,  the  smallest  integer  greater  than  log  (1  C)  -  los  (l  D) 

log  (1-D) 

would  give  required  result.  Similarly,  for  a  given  n  ,  we 
can  easily  determine  the  size  of  the  interval  D  that  would 
be  needed  to  yield  an  expected  coverage  C  . 

Another  question  which  leads  to  a  useful  result  is 
the  following:  given  n+1  points,  chosen  independently  and 
at  random  from  the  unit  interval  (0,1)  (from  which  we  can 
easily  calculate  the  coverage  n(X)  from  (3.1.1)),  what 
would  be  the  expected  increase  in  coverage  resulting  from  the 
random  choice  of  one  more  point,  u  say  ,  from  (0,1)  ? 


nojn-  f  i  s  b  be&ti  ew  tncx)aoup  a irfcf  navrana 


(ci,i;c) 


4.avi«  r  Icfr  *  o  axxe  arW  ©niarta^ot  ^lies*  nao 


:^niwoIXol  eii* 


ei  i  moi  t  x;  I  i»*i  '»j  »  •  roo  n.x  aax>e*ion±  b»&>e<.xe  &d&  &6  bluow 
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Let  AC  denote  the  random  variable  representing  this  increase 

in  coverage  and  let  it  assume  the  values  c  »  In  any  such 

problem  E(AC)  can  be  determined  as  in  the  following  example. 

Suppose  that  we  are  given  that  D  =  1/5  and  that  the  points 

xo  =  0*00  i  x1=0.10  ,  x2~0.35  ,  x^  =  0 ,45  ,  and  x^=0.70 

are  the  ordered  values  of  4  points  which  were  chosen 
independently  and  at  random  from  (0,1)  . 


> - 1 - ( - 1 — | — | - 1 - j - 1 - 1 - 1 

x0  xx  x2  X3  x4 

0  1 
First,  we  calculate  the  distribution  function  of  AC  .  Now, 
since  there  is  a  finite  probability  that  there  is  no  increase 
in  coverage  as  a  result  of  choosing  the  point  u  ,the  proba¬ 
bility  function  of  AC  will  have  a  saltus  point  at  c=0 
of  an  amount  equal  to  Pr( AC =  0) =  0 ,10  +  0.10  =  0.20  . 

Let  0.00  <  c  <  0.05  ,  then 

Pr(AC  <c)  =  o.20+Pr(x1  <  u  <  xL +  c) +  Pr (x2  -  c  <  u  <  x2  ) 

+  Pr(x^  <L  u  <X3+c)+Pr(x^-c<u  <  ) 


+  Pr(x^  <  u  <  x^+c)  +Pr(l  -  c  <  u  <  1) 


' 

. 


ft  .\o:  o  ©*1  >  r  -  a  8 Ji  .  •;  j  "io  i‘.e  v  bei  >t  'io  ■  1  si* 


iO'i 3  t  Lv*ia'i  vi  .  b  :*>  I^fl©  ’flecobni 


Ix 


- 


(I  >  u  >  o  -  C)'t'J+(3+  X  >  u  >  ,  cM  + 
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—  0,2  +  6c  , 

Pr(AC=0.05)  -  0.15  +  0.15  =  0.30  . 

Now,  let  0.05  <  c  <0.10  ,  then 

Pr(AC  <  c)  =  0.80  +  Pr(x^  +  0.05  <  u  <  +  c  )  +  Pr  ( 1  -  c  <  u  <  0.95) 

=  0.70  +  2c  . 

Pr (  AC  =  0.10)  =  0.10  . 

Hence,  the  distribution  function  F(c)  =  Pr(AC  <  c)  of  the 
random  variable  AC  is 


(  0.20 


c  =  0.00 


6c  +  0,20 
F(c)  =  «  0.80 

2c  +  0 . 70 

1.00 


0.00  <  c  <  0.05 
c  =  o .05 


0.05  <  c  <  0.10 
c  =  0.10 


We  readily  obtain,  therefore. 


o.o5 


E(AC)  =  (0.20)  (0.00)  +  6  J  cdc  +  (0.30)(0.05) 


0.00 

0.10 


+  2  J  cdc  +  (O.IOMO.IO) 
o.o5 


=  o.o4 


. 


or  '  *.  ■*  ■  '  ' 
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So,  for  this  particular  example,  we  would  expect  that  the 
random  choice  of  one  more  point  from  the  unit  interval  would 
increase  the  coverage  by  0.04  . 

The  previous  problem  might  lead  one  to  pose  the 
following  question:  given  n+1  points,  which  have  been 
chosen  independently  and  at  random  from  the  interval  (0,1)  , 
on  the  average,  how  many  more  points  would  one  have  to  choose 
from  (0,l)  until  an  increase  in  coverage  is  observed  ?  Let 
N  be  a  random  variable  taking  values  corresponding  to  the 
number  of  points  chosen  until  an  increase  in  coverage  is 
observed.  Returning  to  our  previous  example,  it  is  clear, 
therefore,  that  the  random  variable  N  may  assume  the  values 
1,  2,  3,  ...  with  probabilities  0.8  ,  (0.2)(0,8)  ,  (0.2)2(0.8), 
...  respectively  .  Hence 


00 

E(N)  =  £  k(0.2)k_1(0.8) 

k=l 

00 

=  0.8  £  k(0.2)k-1 

k=l 

=  o.8/(i.o  -  0.2)2  =  5/4  • 


Therefore,  we  would  expect  an  increase  in  coverage  if  two 
more  points  are  chosen. 


io\  r  i  ■  -O'* ^  oq  «  I  :»i  i<  •  «  >  '  ^ 
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ILJL  The  Probability  Distribution  q£  the  Measure  of.  a 

Random  Li-near  Set  for  A  Specified  Number  (Three)  a£ 
Points 

For  the  sake  of  completeness  and  as  a  verification 
of  the  main  result  of  this  thesis,  we  would  now  like  to 
evaluate,  using  methods  quite  different  from  those  of  Votaw, 
the  probability  (and  distribution)  function  of  the  random 
coverage  p(X)  for  the  special  case  n=2  (i.e.  for  3  points). 

From  (3.1.1)  we  have 

(3.2.1)  |i(X)  =  c(x1)  +  c(x2-Xl)  +  c(l-xn)  . 

Consider 

(3.2.2)  E(e^  )  =  J  J  gU-^x^dR  , 

R 


where  g(x-,  ,  x9)  denotes  the  joint  density  function  of  the 

X  X 

order  statistics  and  x2  ,  and  R  denotes  the  region 

over  which  this  density  is  defined.  It  can  be  shown  that  the 
joint  density  function  of  n  order  statistics 
x  ,  x  ,  ...  ,  x  (where  the  density  function  of  each  random 

variable  before  the  ordering  process  is  the  uniform  density 
over  (0,1))  is  nj  for  0  <  xx  <  *2  <  '  ’  *  <  xn  <  1  ‘ 


)o  f  o  ^  1  ) 


■ 


[n  t-i  ((it0) 


5l 


Therefore  (3.2,2)  becomes 

(3.2.3)  E(e^(X)t)  =  2»  J  J  e^c(xl)  +  x2'xl)+ c  ^1-xn^)fc 

R 

where  R  ,  being  the  region  over  which  the  density  function 


g(x1  ,  x2) 


0  <  x!  <  x2  <  1  , 
elsewhere  9 


is  defined,  consists  of  the  triangular  region  in  the  (jc-.  ,  x0 ) 

■1  Lt 

plane  bounded  by  the  lines  x^  =  0  ,  x2  =  1  ,  and  xi  =x2 

However,  due  to  the  possible  overlapping  of  the  intervals, 
as  the  values  of  x^  and  x2  vary  within  this  triangular 
region  R  ,  the  values  of  the  functions  defining  n(X)  , 
namely  c(x1)  ,  c(x2  -  ,  and  c(l-  x2>  (see  (3.1.2)),  also 

change,  so  that  in  order  to  evaluate  (3,2.3)  ,  we  are  forced 
to  integrate  over  complicated  regions  within  R  ,  (it  will 
be  seen  shortly  how  the  number  and  form  of  these  regions 
depend  on  the  value  of  D0) 

Suppose,  now,  that  0  <  D  <  1/3  .  The  form  of  the 

function  c(:x)  dictates  that,  in  order  to  evaluate  (3.2.3)  , 
the  region  R  must  be  decomposed  into  seven  disjoint  regions 
R  (i*l,2,...,7)  ,as  shown  in  the  figure  below. 


dR  , 


. 

- 

. 


J  \  bi<  x,:i  c  fi 


- 
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Notice  that  if  D  =  1/3  ,  the  region  Ry  disappears.  Now, 
since  R  is  the  union  of  disjoint  regions  R^(i = 1,2, . , , , 7)  , 

(3,2,3)  becomes 


(3.2.4)  i  E(eki(X)t)  =  £  J  J  e(o(xl)  +  c(x2-xl)+  c(1-x2})t  dR. 

i=l  Ri  x 


D  x-.  D  .  D  1-D 

x  (x2+D)t  (x  +  2D)t 

=  J  J  e  dx  dXl+J  J  e  1  dx„dx 


0  x. 


0  x1+D 


2  1 


D  1 


1  x2‘D 


p  r*  (x,  +d+i-x„ )t ,  .  ,  r  p  (2D+i-x0)t. 

j  J  e  -1  2  dx2dxl  +  J  J  e  2  dxx 


dx, 


0  1-D 


1-D  D 


■ 
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1  x9  1-2D  x,  +  D  t  > 

+  J  J  (»-x1+l)t  J  J.  e(2D+X2-XlH 

i-D  x2-D  D  X;l  2  1 


1-D  1-D 

+  !  J e 

1-2D  x. 


1-2D  1-D 

( 2D+x  -x-  )t  .  .  ,  p  p  3Dt.  . 

2  1  dx2dx1+  J  j  e  dx2dx1  , 


D  X;l+D 


from  which 


(3*2*5)  =  I2_^0t- l8e^Dt  +  6e.Dt  +  6(1-3D)  e3Dt 


+  6(2D-1)  e2Dt  +  (1-6D+  9D2  )  e3Dt  . 
t 

It  is  easy  to  show  that  lim  E(eU'jl  0)=  1  9  so  that  E(e^^^) 

t  “*  0 

is  indeed  a  moment  generating  function  of  the  distribution  of 
l»i(X)  (for  the  case  0  <  D  <  1/3)  .  Since  the  coefficient  of 
tr/rl  in  the  expansion  of  as  a  power  series  is  the 

r  th  moment  of  |a(X)  ,  it  is  easy  to  show  that 


(3.2*6)  E(n(X))  =  b°  -  3D2  +  3 D  -  1  -  (1-D)J  , 


and 


. 


■ 


•*o  fl0i«AMO  nl 


■ 
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E(h2(X))  =  D4  -  16D3  +  9D2  , 


from  which  we  obtain 


(  3  •  2  •  7  )  Var  ( |i  ( X ) )  = 


-  D6  +  6D5  -  —  D4  +  2D3 

2 


Consider  now  the  function 


(3*2.8)  cp-^x) 


=  /K 


-  to  <  x  <  a  , 

x  >  a  , 


where  K  and  a  are  any  real  constants.  Then 


oo 


(3*2.9)  J*  cp^(x)  ex^  dx  =  Ke 

t 


at 


-GO 


Similarly,  if 


(3.2.10)  <P2(x)  = 


a  K  (  1  -  |  ) 


00  <  x  <  a 


x  >  a  , 


then 


00 


f  cp  (x)  ext  dx  =  K  e_3t 

1  2 

-00  t 


(3.2.11) 
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Now,  an  integral  of  the  form  (3 *2*9)  or  (3.2.11)  can  be 
associated  with  each  term  of  the  expression  (3 .2*5)  (with 
the  exception  of  the  last  term  which  corresponds  to  a  saltus 
point  in  the  probability  function  of  p(X))  .  Therefore, 

(3 *2 *5)  can  be  written  as 


E(e^X>t) 


D 


=  J*  36D(1-  - 


xt 


2D 

f  36  D  (  1- 


D 


x 

2D 


)extdx  +  J*  6D  (1 )extdx 


xt 


“CO 


“CO 


-CO 


+  f  6(l-3D)eXtdx  +  J  6(2D-l)extdx  +  ( 1-6D + 9 D2 )e^Dt 
-co  “°° 

( 3  o  2  o 12 ) 


3D 

(x-D)extdx+  6  f  (-2X+ 3D+l)extdx  +  (  1-6D+ 9D 

2D 


2 


The  term  ( 1-6D +  9D2)e  corresponds  to  a  saltus  or  jump 

point  in  the  probability  function  at  the  point  x=oD  with 
frequency  1-6D+9D2  .  Now,  it  follows  that  the  probability 

function  f2(x)  of  n(X)  is  given  by 


(3.2.33)  f«(x) 


f  0 

6 (x-D ) 

6  ( -2x  +  3D  +  1  ) 
1-6D  +  9D2 


-qo  <  x  <  D  , 

D  <  x  <  2D  , 

2D  <  x  <  3D  , 

x  =  3D  , 

x  >  3D  , 


. 


'•  •  o  >£ 

)c  I 

---.n  an  tm  .  r 


-X 


— ■ 

,  -  *  '  la 


,  l 


provided  0  <  D  <  1/2  •  The  function  given  by  (3*2*13) 

is  continuous  over  the  range  D  <  x  <  3D  with  the  exception 
of  x  =  2D  and  x  =  3D  (x  =  3D  being  a  saltus  point)  • 
Integrating  f^(x)  ,  the  distribution  function  F^(x)  of 
H(X)  is  found  to  be 


r 


(3.2.14)  F2(x)  =  \ 


3 (x-D)2 


-6x2  +  6  (1+  3d)x-9D2  -  12D 


“CD  <  x  <  D 
D  <  x  <  2D 

2D  <  x  <  3D 


v. 


1 


X  >  3D 


which  is  continuous  in  the  range  D  <  x  <  3D  except  at  the 
saltus  point  x  =  3D. 

For  D  =  1/3  ,  the  probability  and  distribution 
functions  may  be  obtained  by  the  same  method  used  above.  As 
mentioned  previously,  for  D =  1/3  ,  the  region  Ry  of  the 
previous  figure  disappears  so  that  we  need  only  integrate 
over  six  regions  in  the  (x^  ,  )  plane.  Equations  (3*2.5) 

(3*2*6)  ,  (3*2.7),  (3.2*13)  ,  and  (3*2*14)  all  hold  for  the 
particular  case  D  =  1/3  *  It  is  interesting  to  note  that 
in  (3*2*13)  ,  the  saltus  point  at  x  =3D  =  1  disappears  and 
the  probability  function  is  continuous  at  this  point.  (This 
saltus  point  was  the  contribution  of  the  region  Ry  to  the 

probability  function  of  n(X)  .)  The  probability  function 
of  KJL ( X )  for  the  case  D  =  1/3  is 


•  r.  '  ■  1  s  '."  f'1 

, 

O  T  l  \ 

O 

('V. 


I  1  '3 


ju  i 

1.  ■  •  -  • 

•  ►  &*  ■  9>- 

r  • 

.  *  3.  ri.  .>1  n 
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(3*2,15)  f2(x)  =*\ 


2(3x-  1  ) 
12( 1  - x ) 


-OD  <  X  <  J  , 

1  ^  ^  2 

3  ^  x  ^  3  , 

4  <  x  <  1  , 

X  >  1 


This  function  is  continuous  over  the  range  D  <  x  <  1 
except  at  the  point  x  =  2D  •  Integrating  (  3*2*15)  ,  we 
obtain 


r 


(3,2.16)  F2(x)  =-< 


|  (3x  -l)2 
-6x2  +  12x  -  5 


-oo  <  x  <  ^-  , 

~<x<~ 

3  ^  x  ^  3  * 

|  <  x  <  1  , 

X  >  1  , 


which  is  continuous  in  the  range  D  <  x  <  1 


In  a  similar  way,  we  can  derive  the  corresponding 
results  for  other  values  of  D  (0<D<1)  .  These  results 

are  summarized  below* 

For  the  case  1/3  <  D  1/2  we  have 


(3.2.17)  E  (M(X) )  =  1  -  (1-D)3 

(3.2.18)  Var  (u(X))  =  -  D6+ 6D5- ^  D^+  2D3  > 


■ 
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( 3  o  2 „ 19 )  f2(x) 


=  4 


\ 


0 

~00 

<  X 

<  D 

6 ( x-D) 

D 

<  X 

<  2D 

6  (  -2x  +  3D  +1) 

2D 

<  X 

<  1 

1-6D  +9D2 

X 

=  1 

0 

X 

>  1 

which  is  continuous  in  the  range  D  <  x  <  1  except  at  the 
points  x  =2D  and  x  =  l  (this  point  being  a  saltus  of  amount 
1-6D+  9D2  ),  and 


(3.2.20)  F2(x)  =-( 


3(x  -D)Z 

•6x2  +6(l  +  3D)x  -9D2-12D 
1 


-00 

<  X 

<  D 

D 

<  X 

<  2D 

2D 

<  X 

<  1 

X 

>  1 

which  is  continuous  except  at  the  saltus  point  x  —  1  • 
The  case  1/2  <  D  <  1  yields 


(3.2.21) 


E(p(X) )  =  1  -  (1-  D)-5  , 


Var  (p(X))  =  -  D6  +  6D5-^D4-12D2  +  4D 


(3.2.22) 


N5  |H 
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r  o 

-CD  <  X  <  D 

f 

(3,2,23) 

6(x  -  D) 

D  <  x  <  1 

f,(x) 

r% 

9 

L 

-  2  +  6D  -  3DZ 

x  =  1 

9 

^  0 

x  >  1 

9 

which  is 

continuous  in  the  range 

D  < 

x  <  1  with  the 

exception 

of  a  saltus  of  amount 

-2  + 

6D  -  3D2  at 

X 

and 

r  0 

1 

8 

A 

<  D  , 

(3.2.24)  F2(x)  =  < 

3(x  -  D)  2 

D  <  x 

<  1  , 

C  1 

X 

>  1  , 

which,  is  continuous  in  the  range  D  <  x  <  1  with  the 
exception  of  the  saltus  at  x  =  l  , 

3  ,  3  The  Probability  Distribution  of  the  Measure  of  £  Random 

Linear  Set  for  an  Unspecified  Number  of  Points  , 

We  now  derive  the  main  result  of  this  thesis.  We 

show,  using  the  methods  of  Votaw,  that  the  probability 

function  f  (x)  of  this  random  variable  p(X)  (which  was 
n 

defined  at  the  beginning  of  this  and  the  first  chapters)  is 
given  by  (1.3)  . 


Cl  1  K  >  X. 


. 


*  ?,  Si  i  a-  rfoltf'w 
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Recall  that  oar  circular  representation  of  the 
problem  was  the  following.  Order  n+1  points  chosen 
independently  and  at  random  on  a  circle  of  unit  circumference, 
as  ( i =  0 , 1 , . . . ,n)  ,  To  each  point,  attach  an  arc  of 

length  D(0  <  D  <  l)  ;  the  point  being  at  the  clockwise  end 
of  the  arc.  Choose  the  beginning  of  any  one  arc,  xq  say, 
as  the  origin  and  measure  distances  in  an  anti-clockwise 
direction  around  the  circumference.  Let  ^ji(X)  be  that 
portion  of  the  circumference  which  is  covered  by  the  n+1 
arcs.  Make  the  transformation 

(3.3.1)  yQ  =  Xo  =  0  *  yi  =  xi  1^2,  ...  ,n)  ,  yn+1  =  1  -  xn  . 

n+1 

It  is  clear  that  y .  >  0  ( i =  1, 2 , . . . ,n+l)  and  Z  y .  =  1  , 

1  i=l 

so  that  we  will  be  concerned  only  in  those  points 
(y1>y2>»««>yn+i^  which  lie  on  that  portion  of  the  hyperplane 

yi+y2+.c0+yn+1=  1  for  which  y±  >  0  , y2  >  0  ,  . .  . ,  ^+i  >  0  . 

Denote  this  region  by  S  .  The  content  of  S  ,  V(S)  say,  is 
given  by  Hallos  result  (see  (203o4-)  )  es 


. 


O'  to  il'i  40  3;,’  r  u  i'  >ci  iwt  .  J*  c*  'i 
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(3.3.2)  V(S)  =  i£±I 

*>! 

If  R  is  any  region  of  S  ,  with  content  V(R)  ,  then  since 
points  (y-^  >y 2 »  9  •  9  >yn+l)  are  e<lua^ly  likely  to  be  at  any 

position  in  S  ,  the  probability  that  a  point  p  belongs  to 
R  is 

v(r)  n  o 

(3.3.3)  Pr ( p  6  R)  =  V(R)  . 

(Note  that  from  (3.3.2)  it  follows  that  the  joint  density 

n+1 

function  of  y^  ,y2 ,  .  *  o  ,yn+^  (y^  —  0  ;  i=l ,  2 , . .  .  ,n+l ;  Z  y^  =  1 ) 
is  just  nl/Jn+1  .) 

Now,  we  have  seen  that  p(X)  is  given  by 

n-1 

*i(X)  =  c(x1)  +  Z  c(xi+1  -  x±)  +  c(l-xn)  , 

i=l 

where  the  function  c(x)  was  defined  by  (3.1.2)  .  The 
range  of  n(X)  is  D  <  ti(X)  <  m  ,  where  m  denotes  the 
minimum  of  1  and  (n+l)D  .  If  m=(n+l)D  ,  there  is  a 
finite  probability  that  the  set  X  will  consist  of  n+1 
disjoint  arcs  of  length  D  , 


and  hence  the  distribution  of 


' 


al 


d&^ael  %o  coin  ^alot»-tb 
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li(X)  will  have  a  saltus  point  at  li(x)  =  (n+l)D  ,  Given  that 
ra  =  (n+l)D  ,  \i(\)  =  (n+l)D  if  and  only  if 

yl  -  D  5  y2  -  D>  0  0  0  »yn+i  —  D  •  Let  R  be  the  subregion  of  S 

consisting  of  all  points  (y  ,y2 , , 0 „ >Yn+±)  which  belong  to 

S  but  also  satisfy  the  condition  that  y^  >  D  ,  >  D, 

. „ . ,  yn+1  >  D  ,  It  is  clear  from  the  diagram  below  that  R 

will  be  the  interior  of  the  convex  polyhedron  whose  vertices 
are 


( 1-nD  ,  D ,  „  *  *  ,  D  ) 
V  (D,l-nD,  ,,OJD) 


Vn+1(D,D, 0 0 ,1-nD) 


(D,D,1-2D) 


<D,1-2D,D) 


(1,0,0) 


1 


( 1-2D, D, D) 


' 

■ 
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We  now  apply  Hall 9 s  result  (namely  equation  (2.3.4)) 
to  evaluate  the  content  of  the  region  R  ,  The  content  of 
this  region  will  be  unaffected  by  translation  by  the  vector 
( -D, -D, . . . , -D)  into  the  polyhedron  with  vertices 


V  * 

1  (l-(n+l)D,  0,  .  .  . , 0) 

(0,l-(n+l)D,...  ,  0) 

a 

o 

o 

vn+l(° ’° ’ ' • • »  l-(n+l)D)  • 

But  these  are  simply  the  points  of  intersection  with  the 
coordinate  axes  of  the  hyperplane  z^+z2+  •  «  •  +  z^+^=l  -(n+l)D 

(z^  =  y^-D;  i  =  1 , 2,  .  .  .  ,n+l)  .  Now,  in  order  to  apply  Hall’s 

result,  we  must  apply  a  uniform  stretch*,  The  substitution 
= z^/(l-(n+l)D)  (i =  1, 2, . . .n+l)  expands  every  dimension 

zi(i= 1,  2,  . . . ,n+l)  by  a  factor  l-(n+l)D  ,  giving  the 

hyper  plane  qi  +  q£  +  °  o  •  +  =  (q^^0;i=l,2,Oo.  ,n+l)  . 

We  have  seen  before  (see  (3.3.2))  that  the  content  of  this 
region  is  given  by  Hall’s  result  (equation  (2.3.4))  as 
Jn+l/nl  .  Contracting  now,  by  the  factor  l-(n+l)D  ,  the 
content  of  the  region  R  is 

(l-(n+l)D)“  . 

“I 


.  ■  •'  *1 


l  I  :  ,  vs'-  3  , 


lojotl  6 
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Hence,  from  (30303)  ,  given  m=(n+l)D  , 


PrMx)  =  (n+l)D)  =  Pr(yx  >  D  ,  y2  >  D ,  yn+1  >  D) 
(3*3.4)  =  (l-(n  +  l)D)n  o 


Hence,  at  the  point  p(X)  =  (n+l)D  ,  the  probability  function 
of  M- ( X )  will  not  be  defined  and  the  distribution  function 
will  have  a  saltus  of  amount  (l-(n+l)D)n  0  Now  consider  the 
distribution  of  p(x)  on  the  range  D  <  p(X)  <  m  (p(X)  takes 
the  value  D  with  probability  zero)  0 


Let  the  random  variable  |i(X)  take  the  values  x 
and  proceed  in  a  manner  analogous  to  the  methods  of  Votaw. 

The  probability  that  X  <  x  <  X  +  AX  (where  X  <  ra  and  AX 
denotes  an  arbitrary  small  positive  increment  in  X  )  can  be 
evaluated  by  determining  volumes  of  certain  regions  contained 
in  S  (the  n-dimensional  region  consisting  of  all  points 
(y^>y2 » • • • >yn+^)  Which  lie  in  that  portion  of  the  hyperplane 


yx+  y,+  ...+  yn+1 


=  1  S 


atisfying  y±  >  0  (i =  1 , 2, 0 c 0 ,n+l)  ) 


Note  that  since  x  <  m  ,  we  can  always  choose  AX  such  that 
X  +  AX  <  m  ,  Consider  the  following  conditions 


(a)  qD  <  x  <  (q+l)D 


( q=l ,2,oo0,M;  M  denotes 
the  minimum  of  n  and  the 

greatest  integer  <  “j  ; 

x  <  m  )  o 


' 


•i 


i 


. 
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Fop  any  particular  value  of  q  ,  decompose  the  set 
^yl  ,y2  9  °  °  0  j  ^  anc*  (y  j  *yj+i >  °  °  °  ?^n+l^  such  ‘that 

(b)  yu  >  D  (u  =1,2,  ,,»,j  ;  j  <q)  , 

(c)  yy  <  D  (v  =  j+l,j  +  2,9S0,n+l)  . 

n-1  n+1 

Now,  we  have  x=c(Xl)  +  E^  c(xi+1~  x± )  +  c(l-xn)  =  E  c(yi)  , 
which, in  view  of  (b)  and  (c)  becomes 


n+1 

(  3  <,  3  »  5)  x  =  3  D  +  E  y  0 

_  .  ,  v 

V  -  J  f  1 


n+1  3  n+1 

Now,  we  have  E  y  =  E  y  +  E  y  =1  ,  which,  with  (3. 3,6) 

i=l  1  u=l  u  v-j+1  v 

becomes 

9 

3 

(d)  E  y  =  1  -  x  +  JD 
u=l 


Let  A.,(x)  denote  the  content  of  the  (j-1)- 

dimensional  region  consisting  of  all  points  (yn *y9 > o • • ,y  . ) 

±  Z  3 

which  satisfy  conditions  (b)  and  (d)  ,  and  let  B  _.(x) 


■ 


• 

■J 


' 
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be  the  content  of  the  (n-j ) -dimensional  region  which  is 
the  intersection  of  the  hyperplane 


n+1 

Z  y  =  x-  jD  with  the 


v=j+l 


(n-j+l) -dimensional  hypercube  given  by  condition  (c) 


The  derivation  of  A.  ,(x)  is  the  same  as  that 

0 


of  (3.3.4) 
of  all  points  (y-^.y, 

y±  >  D ,  y2  >  D 


A.,(x)  is  the  volume  of  the  region  consisting 
3  1 


j  o  o  o ,y  •  )  such  that 
3 


,,,,,y.  >  D  and  which  lie  on  the  (j-l)- 
3 


dimensional  hyperplane  +  ooo+  y  ,  -l-x+jD  0  We  require 

then,  the  content  of  the  (j-l) -dimensional  region  with  vertices 
(each  having  j  components) 


V^(l-x+D,  D,oo0,D) 
v2(d,i-x+d,coo,d) 

o 

o 

o 

V.(D,D,000,l-x  +  D  ) 

3 

Proceeding  as  before,  translation  of  this  region  by  the  vector 
(-D,-D, , o o ,-D)  will  give  a  region  with  vertices 


. 
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V^(l-x,0, 00  o ,0) 

V2^°>1"x>  ««»>°) 

o 

6 

o 

v  (0,0, o  o  o ,1-x)  , 

J 

which  is  the  intersection  of  the  hyperplane  w, +w  +  ,  +w.  =  1-x 

12  j 

(wi  **"  Yi  -  D  j  i  =  1 , 2 ,  o  o  o ,  J  )  with  the  3  -dimensional  hypercube  of 

side  1-x  »  In  order  to  apply  Hall’s  result,  we  let 

3^ =  w^/( l~x) ( i = 1 , 2 , 0 0 o , j )  which  expands  this  region  to  a 

region  with  vertices 


V'(l,0,ooo,0) 

V2  (  0  >  1  9  °  0  °  9  0  ) 


o 


v'(o,o 

3  9 


9 


which  is  the  intersection  of  the  hyperplane  q^  +  +*«*+q^  =  l 

with  the  unit  hypercube*  Using  Hall’s  result,  it  is  clear 
then,  that  the  content  of  this  ( j -1 ) -dimensional  region  is 
Vj/(  j-Dl  •  Thus,  contracting  by  the  same  factor,  the 
required  content  is 


'  -J,..  »'■  f 
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(3'3*6) 


Ad_l(x) 


V  j 

(d-m 


(l-x) 


j-i 


Bn  .(x)  is  the  content  of  the  (n-j ) -dimensional 
region  cut  from  the  hyperplane  y^.  +  -^  +  ^j  +  2  +  ***  +  yn+i  =x-jD 

by  the  hypercube  of  dimension  n~j+l  defined  by 

yj+l  '  "  ’  Jj+2  '">•••  SJ,n+X 


<  D  ,  y^_,_0  <  D  ,  .  .  •  ,  y  <  D  .  In  order  to  apply  Hall’s 


unit  hypercube  by  letting  z^=y^/D  (i=  j  +  1 


result ,  we  expand  the  above  hypercube  of  side  D  to  the 

. . ,n+l)  ,  so 

that,  using  Hall,  we  can  determine  the  volume  of  the 
(n-j ) -dimensional  region  which  is  the  intersection  of  the 

j  with  the  unit 


hyperplane  z j  +  i +  Zj  +  2 +  *  *  *  +  zn+l  =  X/D 
hypercube  2  j  +  l  <  1  *  *  *  *  9  Zn+1  ^  ^  * 


Now  since 


qD  <  x  <  (q-S-l)D 


may  be  written  as 


q-j  <  |  -  j  <  (q+1)  ~  3  , 


we  have  [x/D-j]=q-j  ,  so  that,  from  (2.3*4)  ,  the 
content  of  the  expanded  region  is 


Vn-j+j 

(n-j ) ! 


q-j 

Z  (-D 

r=  0 


r 


n 


j+1 


)  ( 


n-j 


D 
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On  making  the  appropriate  contraction,  we  obtain 


(3.3.7) 


B  .  (x  ) 
n-jv 


Vn-j+1 

(n-j)! 


V(-l)r  (  n-J+1  )[x-D(j+r)]n’3 


Now,  A.  ^(x)  and  B  .(x)  define  contents  of 
J  ^  “  J 

regions  contained  in  orthogonal  subspaces  of  the  n-dimensional 
hyperplane  + ^n+1  = ^  the  region  S  )  . 

Aj_^(x)  is  the  content  of  a  ( j -1 ) -dimensional  region  which 
belongs  to  a  subspace  of  a  vector  space  consisting  of 
(n+1 )  -tuples  (y-i  9y  •••  9y  • ,  0*  •  •  •  *  0  )  while  B  .  (x)  is 

equal  to  the  content  of  an  (n-j ) -dimensional  region  which 
belongs  to  a  subspace  of  a  vector  space  consisting  of 
(ntl)-tuples  ( 0 , . . . ,  0  ,  y^+^,  •  °  •  >yn+i  )  •  These  vector 

spaces  are  orthogonal  compliments  whose  product  space  is  the 
positive  orthant  of  (n+1 ) -dimensional  Euclidean  space*  The 
n-dimensional  region  S  belongs  to  this  product  space  and 

we  wish  to  determine  the  content  of  that  subregion  of  S  , 

T(x)  say,  which  is  the  resulting  product  of  the  regions 
with  contents  given  by  and  Bn_^(x)  °  is > 

T  (x)  is  the  content  of  the  (j-l)  +  (n-j  )  =  (n-1)  -  dimensional 
subregion  of  the  n-dimensional  hyperplane  S  which  consists 
of  all  points  *^2  * ’ * # ,yn+l^  which  satisfy  conditions 

(b)  ,  (c),  and  (d)  (it  is  clear  that,  for  points 


il  a  -  i 


'to  noisdidjji  to  t  9  Lkm  ti  b  oS  ri  *±w  ®w 

5  )  7v  o  >ii$  4  il  il  f  ; 


-  ?o  - 


belonging  to  S  ,  one  of  condition  (d)  or  equation  (3*3.5) 
is  redundant ) .  The  content  of  T(x)  is  therefore, 

Vi<x)  Bn_j(x)  *  which,  from  (3.3.6)  and  (3.3.7)  ,  is 
given  by 

(3.3.8)  T (x)  =  T&TT  (l-x)-*'1  fn-J+1  qSJ(-l)r  (  n’j+1  ) 

(n-j )  l  r=0  r 

•  [x-D( j+r) ]n”J 

Let  us  restate  the  above  conditions  (i.e.  (b)  ,  (c)  ,  (d)  , 


(3.3.5) 

)  as 

(i) 

y1  >  D>  0 • • >  D 

9 

(ii) 

yj+l  <  D'-"’yn+l 

<  D  , 

(iii) 

y  i  +  y 2  +  •  •  •  +  y  .  = 

1  -  x  +  j  D  , 

(iv) 

y  .  .  ,  +  y  + 

3  +  1  3  +  2 

yn+l=x-jD 

We  wish  to  determine  Pr(X<x<X+AX)  ,  given  that  these 

n+1 

conditions  are  satisfied.  Now,  since  (iii)  and  E  y . =  1 

i=l  1 

imply  (iv)  ,  T(x)  is  completely  determined  as  that  region 
of  S  satisfying  conditions  (i)  ,  (ii)  ,  and  (iv)  •  In 
other  words,  consider  the  subregion  S*  of  S  for  which 


Jl<9.X£<  f  d  >  X 
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(i)  and  (ii)  are  true.  Then  T(x)  is  the  intersection  of 
this  n-dimensional  region  S*  with  the  hyperplane,  J(x) 
say,  with  equation  given  by  (iv)  (it  is  clear  that  we  could 
also  have  chosen  J(x)  to  be  the  hyperplane  given  by  (iii)  )  • 
Let  us  consider  an  example.  Suppose  we  choose  3  points  on 
the  circle  of  unit  circumference,  so  that  n+1  =  3  • 

Furthermore,  suppose  j  =  2  •  Then  the  above  conditions 

become 


(i) 

yl  -  D  >  y2  -  D  » 

(ii) 

<  D 

(iii) 

y^  +  =  1  -  x  +  2  D  , 

(iv) 

y3  =  x  -  2  D 

Now,  since  we  want  to  determine  Pr(X  <  x  <  X  +  AX)  ,  we 
must  have  X  <  x  <  X  +  AX  ,  so,  not  only  must  y^  satisfy 
condition  (ii)  ,  but  it  must  also  lie  in  the  region 
X-2D<y^  <  X  +  AX  -  2D  ,  that  is,  between  the  planes  J(X) 
and  J(X  +  AX)  . 


9 


IV 


I  (ifioieaoBjlb~n  B±ri;t 


•olqa&xo  an  mbXmoo  an  J’oJ 
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Hence,  in  order  to  determine  the  probability  that 

X  <  x  <  X  4-  AX  given  that  conditions  (i)  ,  (ii)  ,  (iii)  , 

and  (iv)  are  satisfied,  it  is  clear  from  the  figure  above 

that  we  need  only  evaluate  the  content  of  the  region  ABCE  . 

For  small  AX  ,  this  content  is  given  by  the  length  of  FG 

(which  is  just  the  distance  between  the  planes  J(X)  and 

j(X  +  AX)  measured  along  the  face  of  the  tetrahedron) 

multiplied  by  T(x)  (the  content  of  the  (n-l=l ) -dimin sional 

* 

region  which  is  the  intersection  of  S  with  the  plane 
J(x))  •  It  can  be  shown  that,  in  the  general  case,  the 

length  of  FG  is  AX  */n+l  ,  Therefore,  the  content  of  the 

Vn-j+  1  Jj 

n-dimensional  region  consisting  of  all  points  (y^ >y2  *  *  *  * ,yn+l^ 

satisfying  conditions  (i)  ,  (ii)  ,  (iii)  ,  and  (iv)  and 
the  condition  X  <  x  <  X  +  AX  is,  using  (3»3*8)  ,  given  by 


' 


art  3  ^nacfraoo  arit)  (x)T 
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a-)1-1  pH  P-D-I  )[.-.(j*r>rJ 

Vn-j  +  1  Vj  (  J  -1)  •  (n-j)!  r=o  r  L  J 


Hence,  from  (3«3*3)  ,  dividing  the  above  expression  by  the 
content  of  S  ,  namely  *Jn+\/n\  ,  we  obtain  the  probability 
that  any  point  (y^ ,y  ^ , . . . ,yn+^)  of  S  belongs  to  the 

region  R  consisting  of  those  points  of  S  which  satisfy 
conditions  (i)  ,  (ii)  ,  (iii)  ,  (iv)  ,  and  X  <  x  <  AX  as 


j-1  V<-Dr  <  n-J+1 

r=0  r 


•  [x-D(j+r)] 


in~J 


Now,  if  any  other  of  the  (  n+^  )  combinations  of  the 

j 

j  yfs  had  been  specified  in  (b)  and  the  n-j+1 
complementary  y?s  had  been  specified  in  (c)  ,  the  corres 

ponding  values  of  _^(x)  and  Bn_j(x)  would  be  equal  to 

those  in  (3.3.8)  ,  so  that 


Pr(X  <x  <X  +  4X)  -  Hn!  ^  ^  |n-_^r(j_l). 


(  n+1  )  (  n-J+1  )  (l-x)j-1  [x-D(j+r)] 


n-J 


% 

t  , 

; 

. 

■ 

V, 

■\  ■•  > 

...  •  ( 


^noclfloo 


X+n'Ct  ”  •  * « *«  *•«« 


. — * 


~-x 


Bfi 

!  .  . 


. 
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Hence,  dividing  each  side  of  the  above  equation  by  AX  and 
letting  AX  tend  to  zero,  we  obtain  the  probability  function 
of  |l(X)  as 


f  (x)  =  n» 


n 


q 

E 


q-j 
E 


(-D 


J=1  r=0  (n-j)’.(j-l)'. 


(  n+l  )(  n-J+1  )  (l_x)j-1 


[x-D(j+r)j 


n-J 


for  qD  <  x  <  (q+l)D  ( q=l , 2 , • • • ,M  |  M  denotes  the 

minimum  of  n  and  the 
greatest  integer  less  than 
1/D  ;  x  <  m  )  . 


As  we  have  seen  before,  f^(x)  is  not  defined  at  x=(n+l)D 

if  (n+l)D  <  1  (see  (3*3*4)  )  and  the  distribution  function 
of  |i(X)  is  continuous  with  exception,  in  the  case 
(n+l)D  <  1  ,  of  a  saltus  of  amount  (l-(n+l)D)n  at 
x=(n+l)D  . 


. 


X 


d  ■ 
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